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Abstract 


In  this  report,  we  describe  tools  for  the  creation  and  modification  of 
modular  semi-autonomous  forces  (ModSAF)  terrain  databases  to  support 
the  evaluation  of  a  small  autonomous  robot  in  a  tactical  scenario.  Our  work 
is  motivated  by  the  modeling  and  simulation  needs  of  the  Demo  III  robotics 
program  which  is  developing  a  small  tactical  robot  called  the  experimental 
unmanned  vehicle  (XUV).  The  XUV  is  a  small  wheeled  robot  which  must 
autonomously  navigate  through  its  environment.  The  primary  mission  of 
the  XUV  will  be  to  augment  the  scout  forces,  so  it  must  provide 
reconnaissance,  surveillance,  and  target  acquisition  information  (RSTA)  to 
its  operators.  Modeling  the  XUV  in  a  simulated  environment  is  challenging 
since  existing  terrain  databases  do  not  have  sufficient  resolution  to  examine 
the  mobility  characteristics  of  small  vehicles. 

Our  tools  increase  the  resolution  and  detail  of  existing  terrain  databases  so 
that  the  databases  have  sufficient  detail  to  challenge  the  mobility,  chassis 
dynamics,  and  RSTA  models  of  a  small  unmanned  platform.  To  properly 
model  a  small  vehicle  such  as  the  XUV,  the  terrain  database  in  ModSAF 
needs  to  be  modified.  The  modification  is  done  in  two  phases.  In  the  first 
phase,  the  resolution  of  the  grid  underlying  the  terrain  is  increased  by 
placing  additional  elevation  grid  posts  between  the  existing  posts. 
Elevations  are  assigned  to  the  new  grid  posts  using  mathematical  terrain 
models  such  as  the  variable  resolution  terrain  Model  (Wald  &  Patterson, 
1992).  The  new,  higher  resolution  terrain  directly  affects  the  vehicle 
dynamics  and  the  line-of-sight  algorithms.  The  new  terrain  does  not 
directly  affect  the  ModSAF  route-planning  algorithms.  In  the  second  phase 
of  our  terrain  database  modifications,  the  slopes  on  the  new  terrain  are 
examined.  Regions  that  are  steep  or  inaccessible  to  the  XUV  are  marked  as 
obstacles  in  the  database.  The  route-planning  algorithms  use  these 
“obstacles”  to  avoid  planning  routes  through  regions  that  are  too  steep  for 
the  XUV. 
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MODIFYING  MODSAF  TERRAIN  DATABASES  TO  SUPPORT  THE  EVALUATION 
OF  SMALL  WEAPONS  PLATFORMS  IN  TACTICAL  SCENARIOS 


1.  INTRODUCTION 

In  this  report,  we  describe  tools  for  the  creation  and  modification  of  modular  semi-autonomous 
forces  (ModSAF)  terrain  databases  to  support  the  evaluation  of  a  small  autonomous  robot  in  a 
tactical  scenario.  Our  work  is  motivated  by  the  modeling  and  simulation  needs  of  the  Demo  III 
robotics  program.  Under  the  Demo  III  robotics  program,  the  U.S.  Army  is  developing  a  small, 
survivable,  experimental  unmanned  ground  vehicle  (XUV)  (see  Figure  1)  capable  of  autonomous 
operation  over  rugged  terrain  as  a  part  of  a  mixed  military  force  containing  both  manned  and 
unmanned  vehicles.  The  primary  role  of  the  XUV  will  be  to  augment  the  Army  battalion  and  brigade 
task  forces  scout  platoon.  In  this  scout  mission,  the  XUV  is  expected  to  move  through  the  terrain 
using  proper  military  movement  techniques  and  with  minimal  human  oversight.  Using  its 
reconnaissance,  surveillance,  and  target  acquisition  (RSTA)  package,  the  XUV  can  acquire 
information  about  the  disposition  of  enemy  forces.  We  must  include  “models”  of  XUV  chassis 
dynamics,  XUV  mobility,  and  the  XUV  RSTA  package  to  properly  assess  the  contribution  of  the 
robot  to  the  overall  mission.  It  is  also  important  to  use  terrain  databases  that  have  sufficient  detail  to 
challenge  the  mobility,  chassis  dynamics,  and  RSTA  models. 


Figure  1.  The  XUV. 
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To  properly  model  a  small  vehicle  such  as  the  XUV,  the  terrain  database  in  ModSAF  needs 
to  be  modified.  The  modification  is  done  in  two  phases.  In  the  first  phase,  the  resolution  of  the 
grid  underlying  the  terrain  is  increased  by  placing  additional  elevation  grid  posts  between  the 
existing  posts.  Elevations  are  assigned  to  the  new  grid  posts  using  mathematical  terrain  models 
such  as  the  variable  resolution  terrain  model  (Wald  &  Patterson  1992).  The  new,  higher  resolution, 
terrain  directly  affects  the  vehicle  dynamics  and  the  line-of-sight  (LOS)  algorithms.  The  new 
terrain  does  not  directly  affect  the  ModSAF  route-planning  algorithms.  In  the  second  phase  of  our 
terrain  database  modifications,  the  slopes  on  the  new  terrain  are  examined.  Regions  that  are  steep 
or  inaccessible  to  the  XUV  are  marked  as  obstacles  in  the  database.  The  route-planning  algorithms 
use  these  “obstacles”  to  avoid  planning  routes  through  regions  that  are  too  steep  for  the  XUV. 

Many  ModSAF  terrain  databases,  which  the  developers  refer  to  as  compact  terrain 
databases  (CTDBs),  have  elevation  posts  spaced  uniformly  125  meters  apart.  At  this  resolution, 
it  is  difficult  to  examine  mobility  issues.  As  an  example,  consider  the  effect  of  terrain  slope  on 
cross-country  travel  at  Fort  Hood,  Texas.  Figure  2  gives  the  change  in  elevation  from  one 
elevation  post  to  the  next  required  to  produce  the  desired  slope  across  the  elevation  grid  square. 
The  change  in  elevation  required  to  produce  a  30  slope  on  a  125-meter  grid  square  (72.2  m)  is 
more  likely  to  occur  in  mountainous  regions  where  wheeled  vehicles  are  not  likely  to  go.  The 
125-meter  resolution  database  representing  Fort  Hood  is  almost  flat  in  most  spots,  presenting  no 
significant  slopes  to  climb.  However,  Figure  3  illustrates  that  even  within  a  single  grid  square, 
there  may  be  significant  slopes  to  climb  such  as  those  associated  with  eroded  areas,  ditches,  or 
culverts.  With  the  length  of  the  vehicle  being  used  to  estimate  the  dimensions  of  the  terrain  area 
shown  in  the  picture,  the  photograph  shows  a  1 00-  by  1 00-meter  area  of  the  terrain.  This  area 
would  be  represented  in  most  terrain  databases  as  a  single  flat  square. 


0  20  40  60  80 
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Figure  2.  Change  in  Elevation  as  a  Function  of  Climb  Angle. 
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Figure  3.  A  100-meter  Patch  of  Terrain  at  Fort  Hood.  Texas. 


There  is  information  in  a  CTDB  database  concerning  trafficability.  Each  grid  square  has  a  soil 
type  associated  with  it.  By  using  a  program  such  as  the  North  Atlantic  Treaty  Organization 
(NATO)  reference  mobility  model  (NRMM),  (Ahlvin  1992),  users  can  adjust  the  maximum  speed 
of  their  vehicle  for  each  soil  type.  In  addition,  there  is  an  abstract  layer  in  CTDB  databases  that 
consists  of  polygons  representing  buildings,  lakes,  tree  canopies,  impassable  regions,  and  other 
features  on  top  of  the  actual  terrain  skin.  These  polygons  interact  with  the  route-planning 
algorithms  by  providing  potential  obstacles.  They  also  can  interact  with  the  visibility  algorithms 
by  blocking  the  LOS  between  points.  As  far  as  mobility  is  concerned,  there  are  only  two  possible 
responses  to  each  type  of  abstract  region:  a  class  of  vehicle  may  avoid  all  polygons  of  a  given  type 
or  may  ignore  all  polygons  of  that  type. 

This  behavior  is  a  little  simplistic  for  vehicles  that  travel  through,  instead  of  around,  the  tree 
canopies.  Canopy  polygons  do  not  affect  the  speed  or  the  route  of  the  vehicles  that  travel  through 
them.  In  addition,  all  canopy  polygons  on  a  given  terrain  database  have  the  same  effect  on  vehicle 
mobility.  An  important  parameter  of  the  canopy  that  is  not  considered  for  mobility  evaluation  is 
tree  density  (density  is  considered  in  the  visibility  analysis).  Cultivated  tree  farms  and  orchards, 
with  uniformly  spaced  trees  and  little  underbrush,  should  be  easier  to  move  through  than  natural 
wooded  areas  with  variably  spaced  trees  and  lots  of  underbrush. 

In  evaluating  the  potential  contribution  of  the  XUV  to  the  scout  mission,  it  is  important 
that  we  consider  the  impact  of  the  XUV’s  mobility  on  the  mission.  This  is  done  by  first 
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increasing  the  resolution  of  the  terrain  database.  The  information  in  the  high-resolution  terrain 
database  is  used  to  adjust  the  soil  types  and  to  create  abstract  regions  to  decelerate  and  redirect 
the  XUV.  This  report  presents  only  the  methods  used  to  modify  the  ModSAF  terrain  databases 
to  make  them  more  useful  for  the  evaluation  of  small  vehicles.  Actual  evaluations  of  the  XUV 
will  be  presented  in  future  reports; 


2.  INCREASING  THE  RESOLUTION  OF  THE  ELEVATION  GRID 

Increasing  the  resolution  of  the  elevation  grid  is,  in  theory,  a  relatively  easy  process.  We 
simply  need  to  find  a  higher  resolution  terrain  database  for  our  battlefield  and  use  it  to  construct 
the  new  elevation  grid.  However,  it  is  not  always  possible  to  find  such  a  database.  Standard 
National  Imagery  and  Mapping  Agency  digital  terrain  elevation  data  (DETED)  databases 
covering  most  of  the  world  use  100-meter  grid  posts.  As  a  substitute  for  measured  elevation 
data,  we  can  artificially  increase  the  resolution  of  the  database  by  fitting  the  original  elevation 
posts  with  a  mathematical  function  that  fits  the  measured  elevation  posts  to  the  desired  accuracy 
and  gives  statistically  realistic  surface  variations  between  the  measured  posts.  One  such  function 
is  variable  resolution  terrain  (VRT)  model. 


The  VRT  model  is  a  continuous,  differentiable  surface  generated  by  summing  several 
simpler  surfaces  referred  to  as  hills.  The  equation  of  a  hill  function  defined  on  the  two- 
dimensional  space  of  real  numbers,  9l2,  is  written  as 


h(x,y)=  ae 


f(ll(x,y),(^)lld) 


(1) 


Here,  a  e  %(%,  77)  e  Si2  and  || (x,y),(£,  tj)[|  is  a  metric  on  9t2.  The  most  familiar  metric  is  the 

Euclidian  distance  between  points.  The  Euclidian  metric  gives  relatively  smooth  hills,  so  other 
metrics  are  often  used  to  produce  hills  with  sharper  peaks. 


By  varying  the  metric,  the  function  f,  and  the  parameter  a,  it  is  possible  to  generate  hills  of 
any  size  and  shape.  To  create  a  generic  VRT  surface,  hills  of  various  sizes  and  shapes  are 
combined  using  the  principle  of  self-similarity.  Self-similar  objects  are  invariant  with  respect  to 
scale  so  that  a  portion  of  the  object,  viewed  at  the  proper  magnification  resembles  the  whole 
object.  In  a  generic  VRT  surface,  the  distribution  of  hills  is  statistically  self- similar.  More  details 
of  the  VRT  model  are  given  in  the  paper  by  Wald  and  Patterson  (1992)  and  the  later  papers  by 
Wald  (1994, 1995).  These  later  papers  discuss  methods  used  to  fit  existing  terrain  databases 
with  VRT  surfaces.  The  same  software  tools  we  describe  in  this  report  can  also  create  new 
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terrain  databases  for  ModSAF  with  specific  characteristics  designed  to  test  movement  algorithms 
developed  to  model  the  movement  of  small  vehicles  within  a  ModSAF  exercise. 

The  software  we  have  developed  is  based  on  software  and  documentation  available  in  the 
ModSAF  developers’  kit  (Braudaway  et  al.  1996).  CTDBs  use  compression  methods  to  store 
terrain  databases  in  much  less  space  than  other  terrain  database  formats  require.  In  gridded 
databases,  the  information  about  the  terrain  surface  is  stored  independently  from  detailed 
information  about  the  feature  layer  on  top  of  the  surface.  For  each  elevation  post,  the  elevation, 
two  soil  types,  and  flags  indicating  the  presence  or  absence  of  features  such  as  buildings,  trees, 
roads,  or  canopies  are  stored  in  32-bit  words.  The  elevation  grid  is  broken  into  patches  (a  square 
region  typically  four  posts  long  by  four  posts  wide).  Each  patch  contains  detailed  feature 
information,  such  as  the  location  and  size  of  buildings  or  the  location  and  width  of  roads. 

Features  that  pass  through  more  than  one  patch  must  be  subdivided  so  that  a  portion  of  the 
feature  is  stored  on  each  patch. 

Much  of  our  program  add_vrtc  (given  in  Appendix  A)  is  based  on  the  program  recompile. c 
in  the  CTDB  library  in  the  ModSAF  developers’  kit.  Recompilex  converts  older  CTDB 
databases  to  the  current  format.  It  can  also  add  features  such  as  additional  roads,  bridges,  or 
canopies  to  CTDB  databases.  We  have  extended  this  code  so  that  we  can  change  the  resolution 
of  the  CTDB  databases.  With  our  new  program,  we  can  change  the  soil  types  of  the  elevation 
posts  to  affect  the  movement  of  the  vehicles  on  that  post,  simulating  rough  or  rocky  regions  on 
the  terrain.  Impassable  regions  can  be  added  to  the  databases,  especially  inside  the  existing 
canopies,  thus  forcing  vehicles  to  modify  their  plans  to  avoid  the  impassable  regions. 

Increasing  the  resolution  of  a  database  is  relatively  easy.  The  C  routines 
get_original_elevations,  look _for Jeatures,  and  look _for_canopy_and_lakes,  store  information 
from  the  original  database  so  that  it  can  be  used  to  construct  the  new  database.  A  major 
difference  between  our  code  and  recompilex  is  that  we  must  recompute  the  feature  flags  for  each 
post  in  the  new  database.  For  example,  suppose  a  road  intersects  a  post  on  the  original  125- 
meter  database.  At  25  meters’  resolution,  this  post  is  covered  by  a  5  by  5  square  of  25-meter 
posts.  It  is  unlikely  that  the  road  will  intersect  all  25  of  these  smaller  posts.  These  flags  are  used 
to  streamline  search  routines  for  the  LOS  and  planning  algorithms;  it  is  not  advisable  to  simply 
pass  the  flags  from  the  old  posts  to  the  new  ones. 

Our  fitting  routine  is  contained  in  the  C  routine  make_vrt_hills.  This  routine  uses  the 
fitting  method  outlined  by  Wald  (1994).  In  his  method,  the  VRT  surface  is  determined 
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iteratively.  The  highest  elevation  post  on  the  database  determines  the  height  and  location  of  the 
first  hill.  The  parameters  for  the  first  hill  are  chosen  to  minimize  the  differences  between  the 
elevation  grid  and  this  hill.  Next,  we  compute  the  differences  between  the  elevation  posts  and  the 
first  hill.  The  largest  of  these  differences  determines  the  height  and  location  of  the  second  hill. 

The  process  continues  until  the  difference  between  the  each  of  the  elevation  posts  and  the  height 
of  the  VRT  surface  at  that  point  is  small  enough. 

To  produce  a  new  terrain  database,  we  must  write  the  new  elevation  grid  and  rewrite  all  the 
feature  information.  The  routine  write_elevation _posts  transfers  the  refined  elevation  grid  to  the 
new  database.  The  C  functions,  collect  Jrees,  collect jcanopies,  collect Jbuildings,  and 
collect  linears,  allow  us  to  collect  the  features  from  the  original  database.  These  routines  were 
modified  from  their  original  form  in  recompiles,  so  that  features  could  be  properly  subdivided 
for  the  new  patch  size.  Figure  4 A  shows  a  contour  map  of  a  10-kilometer  by  10-kilometer 
section  of  Fort  Hood,  sampled  at  a  125-meter  resolution.  The  contours  are  equally  spaced,  2 
meters  apart.  Figure  4B  shows  the  same  section  enhanced  with  VRT  hills,  sampled  at  50  meters’ 
resolution.  The  enhanced  database  has  a  number  of  small  hills  and  valleys  that  are  not  present  in 
the  original  database.  The  heights  of  the  original  elevation  posts  have  been  preserved. 


A  B 


Figure  4.  A  10-kilometer  bv  10-kilometer  Section  of  the  Fort  Hood  Database  (2 -meter  contours). 
(A  shows  the  original  125-meter  database;  B  shows  the  VRT-enhanced  50-meter 
database.) 
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Using  this  program,  it  is  possible  to  create  very  high-resolution  databases.  However,  there 
is  a  trade-off  between  database  resolution  and  usability.  Tests  run  at  Fort  Knox  (Nida,  1998) 
indicate  that  it  is  not  practical  to  use  extremely  high  resolution  databases  in  large  ModSAF 
exercises  involving  hundreds  of  simulated  entities.  The  amount  of  time  ModSAF  spends 
processing  terrain  information  increases  to  the  point  that  operators  must  wait  a  long  time  for 
tactical  display  screens  to  be  redrawn.  Also,  the  linear  networks  that  describe  the  roads,  rivers, 
and  other  linear  features  become  unwieldy  since  there  must  be  at  least  two  points  on  each  patch 
that  intersects  them.  The  planning  algorithms  limit  the  number  of  path  points  that  can  be  stored 
for  a  given  plan.  This  number  can  be  increased  but  not  to  the  point  that  it  can  do  much  good. 

At  this  time,  there  are  four  solutions  to  this  problem.  First,  we  can  limit  the  dimensions  of 
the  battlefield.  The  Fort  Knox  tests  involved  databases  that  represented  a  50-kilometer  by  50- 
kilometer  section  of  Fort  Hood,  Texas.  In  its  original  125-meter  resolution,  the  database  was 
represented  by  a  400  by  400  elevation  grid  and  10,000  500-meter  patches  containing  feature  data. 
At  50  meters’  resolution,  the  elevation  grid  was  a  1000  by  1000  matrix  and  the  feature  data  were 
contained  in  62,500  200-meter  patches.  By  reducing  the  dimensions  of  the  battlefield,  the 
amount  of  data  becomes  easier  to  handle.  This  solution  is  well  suited  for  studies  involving  small 
vignettes  in  which  the  entities  are  confined  to  a  small  portion  of  the  battlefield. 

The  second  possible  improvement  is  to  increase  the  dimensions  of  the  patch.  This  option 
was  discarded  since  it  affects  the  portability  of  databases  that  we  create.  A  potential  user  would 
be  forced  to  reset  a  parameter  in  the  ModSAF  code  to  use  the  databases.  Also,  the  dimensions  of 
a  patch  were  chosen  by  the  ModSAF  developers  to  optimize  memory  usage.  Changing  the 
dimensions  may  decelerate  rather  than  accelerate  the  LOS  and  route-planning  algorithms. 

The  third  option  is  to  use  ModSAF  micro-terrain  to  add  more  detail  to  the  terrain  surface. 
Micro-terrain  is  a  layer  of  triangular  facets  on  top  of  the  terrain  surface.  It  is  used  to  describe 
small  terrain  features,  such  as  single  boulders  or  berms.  It  is  also  used  to  describe  multilevel 
surfaces  such  as  underpasses  and  overhangs.  There  are  limits  on  the  micro-terrain  approach. 

First,  any  changes  should  stay  within  a  single  patch  so  the  added  hills  must  effectively  vanish 
outside  this  patch.  There  is  also  a  limit  on  the  number  of  triangles  that  can  be  added  to  a 
database.  The  fourth  option  is  to  create  databases  based  on  triangularized  irregular  networks 
(TINs)  rather  than  uniformly  spaced  grids.  The  triangles  that  comprise  the  TIN  are  not  uniform 
in  size.  Large  triangles  are  used  in  relatively  flat,  featureless  regions  of  the  terrain.  Small  triangles 
are  used  in  rough  terrain.  The  non-uniform  size  allows  posts  to  be  concentrated  where  greater 
detail  is  required.  The  TIN  representation  of  terrain  is  actually  an  extension  of  the  method  used 
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to  add  micro-terrain.  In  TIN  databases,  there  is  no  elevation  grid  and  the  entire  surface  is 
represented  by  triangular  facets.  We  are  still  examining  this  option. 

3.  ADDING  ABSTRACT  REGIONS  AND  ADJUSTING  SOIL  TYPES 

Refining  the  elevation  grid  in  a  CTDB  database  increases  the  probability  that  a  ModSAF 
vehicle  will  encounter  a  significantly  steep  region  of  the  terrain  as  it  travels  cross  country. 
Unfortunately,  the  elevation  grid  is  not  used  by  the  route-planning  algorithms  so  the  vehicle 
cannot  avoid  a  steep  spot  that  is  represented  only  by  the  elevation  grid.  In  ModSAF,  the 
planning  routines  use  the  abstract  layer  of  the  terrain  database  to  plan  the  route.  This  layer 
consists  of  polygons  representing  features  such  as  tree  canopies  and  lakes,  and  line  segments 
representing  features  such  as  railroads,  power  lines,  and  fences.  By  using  the  simplified  terrain 
representation  in  the  abstract  layer,  the  route-planning  algorithms  efficiently  plan  routes  that 
avoid  these  regions. 

We  have  to  add  polygons  to  the  abstract  layer  of  the  CTDB  database  to  force  the  vehicle  to 
“see”  the  steep  regions.  The  routine  add_steep_regions,  shown  in  Appendix  A,  is  used  to  add  steep 
regions  to  the  terrain  database.  This  routine  uses  the  dimensions  of  the  VRT  hills  to  assign  a 
“steepness”  value  to  the  entire  hill.  Steep  VRT  hills  are  enclosed  in  abstract  soil  regions.  These 
abstract  soil  regions  are  simply  polygons  with  an  additional  parameter  that  specifies  a  soil  type  for 
the  entire  region.  Typically,  these  regions  are  used  to  designate  lakes  in  the  database.  By  specifying 
another  soil  type,  we  can  use  these  abstract  soil  regions  to  control  the  movement  of  vehicles  on  dry 
land.  By  altering  the  vehicle  parameter  file,  it  is  possible  to  force  specific  types  of  vehicles  to  avoid 
these  regions.  An  example  of  a  vehicle  parameter  file  is  shown  in  Appendix  B.  In  the  terrain 
analysis  section  of  the  entity  parameter  file,  the  avoidance  mask  includes  the  TERRAIN_WATER 
flag  which  directs  the  vehicle  to  avoid  soil  regions.  The  specific  types  of  soil  regions  are  given  in  the 
avoid  soils  subsection.  In  this  section,  the  soil  type  is  referenced  by  a  binary  code.  This  32-bit 
integer  contains  mobility  parameters  and  the  Simnet  and  close  combat  tactical  trainer  (CCTT)  indices 
for  the  soil  type.  Table  1  shows  a  map  of  the  bits  for  the  binary  soil  code. 


Table  1.  Bit  Definition  for  the  Binary  Soil  Code 


Bits  31-13 

12 

11 

10 

9 

□□□□□ 

BBDQ 

Road 

(R> 

E551 

IfRBgj 

ESI 

WEm 

Unused 

Mobility  parameters 

CCTT  soil 
index 

Simnet  soil 
index 
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Table  2  gives  the  indices,  description,  and  binary  code  for  the  16  most  common  soil  types. 
Generally,  CTDB  databases  use  the  Simnet  soil  types,  although  it  is  possible  to  define  additional 
soils. 


Table  2.  ModSAF’s  Soil  Types 


Simnet 

index 

CCTT 

index 

Soil  type  description 

Mobility 

parameters 

Binary 

code 

R 

S 

N 

W 

0 

0 

Undefined 

0 

1 

20 

Pavement,  concrete,  asphalt 

* 

4417 

2 

17 

Dry  loose  surface  road 

* 

4370 

3 

5 

Soft 

83 

4 

28 

Water  60-inch  depth 

* 

* 

1988 

5 

26 

Water  16-inch  depth 

* 

* 

2981 

6 

4 

Very  soft 

* 

2118 

7 

18 

Wet  loose  surface  road 

* 

* 

6439 

8 

15 

Very  hard  (slippery) 

* 

2296 

9 

19 

Swamps,  bogs — slow  go 

* 

2361 

10 

7 

Hard 

2170 

11 

7 

Hard 

123 

12 

7 

Hard 

124 

13 

22 

Brush  land 

* 

2413 

14 

9 

Very  hard 

* 

1182 

15 

19 

Swamps  bogs — no  go 

* 

1343 

To  investigate  the  performance  of  the  XUV  in  a  scout  mission,  it  may  be  necessary  to  allow 
the  vehicles  to  move  through  the  tree  canopies.  Tree  canopies  are  polygonal  regions  in  the  abstract 
layer;  they  can  be  treated  like  obstacles  by  a  class  of  vehicles  or  ignored.  Treating  the  canopies  as 
obstacles  is  realistic  for  large  vehicles.  Ignoring  the  canopies  so  that  vehicles  can  move  through  the 
areas  of  the  battlefield  covered  by  canopies  is  not  realistic.  It  assumes  that  vehicles  can  negotiate 
any  region  of  the  canopy  at  the  same  rate  of  speed.  The  routine  look _for_canopy_and._lakes  can 
be  used  to  adjust  the  soil  parameters  of  the  underlying  posts  within  the  canopies.  By  changing  the 
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mobility  parameters  in  the  vehicle  parameter  file,  vehicles  will  decelerate  inside  the  canopies. 
However,  the  vehicles  will  still  drive  in  straight  lines.  In  the  mobility  parameters  section  of  the 
vehicle  parameter  file  given  in  Appendix  B,  the  vehicle  decelerates  on  water  (soils  4  and  5)  and  on 
soils  9  through  15.  In  this  section,  the  soils  are  referenced  by  their  Simnet  soil  index  which  is 
shown  in  Table  2.  By  adding  obstacles  to  the  canopies,  the  vehicles  will  use  realistic  routes 
through  the  canopies.  The  routine  add_subcanopies  alters  the  canopies  by  adding  artificially 
generated  impassable  regions  to  the  canopies.  Just  like  the  steep  regions,  the  impassable  regions 
are  abstract  soil  regions — only  the  soil  parameter  is  different.  An  example  of  an  altered  canopy 
regions  is  shown  in  Figure  5.  The  original  canopies  are  represented  by  the  light  gray  hatched 
polygons  in  Figure  5a  and  5b.  The  sub-canopy  regions  shown  as  dark  gray  polygons  in  Figure  5b. 
In  this  particular  example,  the  sub-canopies  are  relatively  large  but  they  do  not  need  to  be.  CTDB 
databases  do  not  have  limits  on  the  size  of  the  polygons;  however,  there  are  limits  on  the  total 
number  of  vertices. 

As  an  alternative  to  the  sub-canopies,  we  can  replace  the  canopies  with  distribution  of  trees 
and  tree  lines.  This  approach  is  ideal  for  small  battlefields  since  it  affects  both  the  LOS  and  the 
mobility  algorithms.  On  large  high-resolution  battlefields,  it  adds  a  significant  number  of  features 
to  the  patch  data. 
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Figure  5.  Adding  Subcanopies  to  a  CTDB  Database.  (Figure  5a  shows  the  original  ModSAF 
canopies;  Figure  5b  shows  subcanopies  in  dark  gray.) 
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4.  CONCLUSIONS 


We  have  designed  tools  to  modify  CTDB  terrain  databases  to  support  the  evaluation  of  a 
small  autonomous  robot  in  a  tactical  scenario.  The  first  tool  increases  the  resolution  of  the 
elevation  grid.  This  modification  has  the  most  effect  on  the  vehicle  dynamics  model  and  the  LOS 
algorithms.  The  second  set  of  tools  adds  polygonal  regions  to  the  abstract  layer  of  the  terrain 
database.  These  modifications  affect  the  planning  algorithms. 

We  have  begun  to  use  these  databases  in  our  examination  of  the  XUV  in  tactical  scenarios. 
Tests  were  conducted  at  Fort  Knox  during  the  summer  of  1998  to  determine  the  feasibility  of 
using  high  resolution  databases  in  large  ModSAF  exercises.  Although  the  results  were  somewhat 
disappointing,  they  have  given  us  some  useful  information  about  the  use  of  high  resolution 
databases.  By  using  ModSAF  micro-terrain  or  a  TIN  database,  we  may  be  able  to  increase  the 
resolution  of  the  elevation  grid  substantially  without  the  problems  of  the  gridded  approach. 

By  using  steep  regions  and  the  sub-canopy  regions,  we  can  generate  more  realistic  paths  for 
the  vehicles.  However,  these  regions  are  based  on  artificial  modifications  of  the  terrain  database. 
It  is  important  that  we  consider  several  statistical  variations  of  the  same  battlefield  in  a  study  to 
avoid  results  that  hinge  on  the  exact  location  or  shape  of  these  regions. 


INTENTIONALLY  LEFT  BLANK 
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APPENDIX  A 
THE  ADD  VRT.C  CODE 
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THE  ADD  VRT.C  CODE 


A.1  Globals 

Globals 
#ifndef  lint 

static  char  rcsid  0  =  "$RCSfile$  $Revision$  $State$"; 

#endif 

#define  DEBUG 
#define  READMES 
#include  <sys/types.h> 

#include  <sys/stat.h> 

#include  <sys/mman.h> 

#include  <fcntl.h> 

#include  <libcmdline.h> 

#include  <math.h> 

#ifndef  sgi 

#  ifhdef  MAP_AUTOGROW 

#  define  MAP JUJTOGROW  0 
#  endif 
#endif 

#define  TMPFILE  ,,recompiIe.tmp,, 

#define  PHY  SICAL_FE  ATURE_SIZE  40*1024*1024 
#define  MAX_NUM_LINEAR_MODELS  400000 
#define  MAX_NUM_AGGREGATE_MODELS  50000 
#define  DEFAULT  TRUNKJRADIUS  0.5  /*  meters  */ 

#define  DEFAULT_FULLNESS  0.95 
#defme  NOFEATURES  -1 
#include  <stdio.h> 

#include  <ermo.h> 

#include  <stdstring.h>  /* common/include/global*/ 

#include  <libctdb.h> 

#include  <ct_feat.h> 

#include  "ct_cmplr.h" 

#include  "ct_reform.h,, 

#include  <ct_post.h> 

#include  <libreader.h> 

#include  <stdalloc.h> 

#include  <stdext.h> 

#include  <time.h> 

#include  <string.h> 

#include  <libcoordinates.h> 

#include  <libgcs.h> 

#include  <curses.h> 

#define  my_factor  1 
int32  myjpatch_x,myj)atch_y; 
struct  MY  HDR  {  float64  incrD,inv_incr,patch_incrD 
float64  max_z,min_z,  patch_inv_incr7post_increment; 
int32  incr,patch_incr,max_x,min_x,max_y,min_y,max_x_post, 
int32  min_x_post,  min_y_post,max_y_post; 
int32  pages_wide,pages_high,patches__wide,patches_high;}  my_hdr; 
#define  INPUT_MEMORY  0 
#define  START_CELLS  50 
static  CTDB_CMP  *ctdb_in,  ctdb_out; 


The  boldface  type  indicates  include  files  that  can  be  found  in  the  ModSAF 3.0 
source  code . 
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static  int32  output_format; 
static  int32  input_format; 
static  int32  feature_offset  =  0; 

static  int32  mf;static  CTDB_COMPILER_QUAD  quad_root; 
static  int32  num_corrections  =  0; 
static  READER_UNION  *corrections; 
static  int32  output_db_patches_wide  =  0; 
static  int32  output_db__patchesJiigh  =  0; 
static  int32  output_db_Jncr  =  0; 
static  int32  output_db_patch_incr  =  0; 
static  int32  offset_x  =  0,  offset_y  =  0; 
static  int32  offset_x_meters  =  0,  offset_y_meters  =  0; 
static  int32  min_patch_x,  max_patch_x,  min_patchry,  max_patch _y; 
static  int32  recompile_buildings  =  0; 
static  int32  recompilejinears  =  0; 
static  int32  recompile_trees  =  0; 

static  float64  indb_south_lat,  indbjiorthjat,  indb_westJong,  indb__eastJong; 
static  int32  gc$_mode; 
static  COORD_TCC_PTR  tcc; 
static  int32  cur_cell; 
static  int32  water_fix; 

static  CTDB  JDBASE_TYPE  dbasejype  =  CTDBJHYBRID; 
static  int32  te_offset=0; 
static  int32  total_tes=0; 
static  int32  trans_end  =  0; 
static  FILE  *fptr,  *color_flle; 

static  TRANSITION J>ATCHJNFO  *trans_patches  =  NULL; 
static  void  look_forjeatures(); 
static  void  add_canopies(); 
static  int32  encode_physicaI_features(); 
static  void  collect_nodes_and_edges(); 
static  void  encode_abstract_features(); 
static  void  complete_micro_poly(); 
static  void  collect_microterrain(); 
static  void  collect_buildings(); 
static  void  collect_trees(); 
static  void  collect_canopies(); 
static  void  collect_Jinear(); 
static  int32  *transform_posts(); 
static  void  fit_vrt(); 
static  double  vrt(); 
static  double  singlejhill(); 

#define  CMP_MAX_REF  1000 
#defme  CTDB_INVALID_REF  -1 
typedef  struct  cmp_ref 
{ 

int32  user_id; 
int32  index; 

}  CMP_REF; 

struct  {  CMD_STRING_OPTION  input_ctdb; 

CMD_STRING_OPTION  output_base_name; }  options  = 

{ 

{"Input  Ctdb",  "Absolute  pathname  to  the  input  CTDB"  , 

NULL ,  CMD  STRING ,  "input_ctdb"  ,  NULL  ,  "knox-03 1 1  .ctb" , 

{"Output  Name  Base",  "Output  database  base  name"  ,  NULL  ,CMD_STRING ,  "output_base_name" , 
NULL, NULL,} 


In  this  section  of  code,  the  bold  print  specifies  a  function 
prototype .  The  code  for  these  functions  is  found  later  in  this 
appendix. 
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}; 


struct  PARAMETERS  {double  cx,cy,hgt,angle,wgtx,wgty,power,wl,w2;} 
hill[ 1 000]  ,patch Jiills[  1 00]  [  1 00]  [10]; 
int  nhills,  savedj3agesjhigh,saved_pages_vvide; 
float64  hgt_max,original[500][500]  =  {  0.0}; 
float64  temp[500][500],  mat_step  =  10.0; 
short  intmy_soil[500][500]  =  {0,canopy[5500][5500]  =  {0}; 
short  int  hills_per__patch[100][100]; 

A.2  Main 

int  main  ( int  argc,  argv_t  argv) 

{ 

int32  result; 

CTDB_FILE_HEADER_CMP  hdr; 

char  **correction_files,  first_db_name[CTDB_NAME_LENGTH],  *input_db_name; 

int32  num_correction_files,  header_assigned  =  FALSE, i; 

int32  *elev_ptr  =  ctdb_out.elevations,*soil__ptr  =  ctdb_out.soils; 

extern  int  atoi(); 

char  dbname,s; 

int32  leftover_argc; 

argv_t  leftover_argv; 

char  *rest_of_str,  *cpy_rest_of_str,  *current_file; 
int32  err,  chars__processed; 
char  sp  =  ’\040'; 

int32  *cells  =  NULL,ncells,celI Jdx; 

float64  cell_south_lat,  cell_north_lat,  cell_west_long,  cell_east_long; 

float64  cell_min_gcs_x,  cell_max _gcs_x,  cell_min_gcs_y,  cell_max _gcs_y; 

float64  cell_min_utm_x,  cell_max_utm_x,  cell_min_utm_y,  cell_max_utm_y; 

float64  cell_origin_x,  cell_origin_y; 

int32  *patch_water_state  =  NULL; 

int  ix,iy  tloc; 

float64  loc_x,loc_y; 

uint32  p; 

int32  norm_x,  norm__y; 
bzero((void*)&quad_root,  sizeof(quad_root)); 

cmd_process_options(argc,  argv,  &leftover_argc,  &Ieftover_argv, 
(CMD_OPTION  *)&options,  sizeof(options),  TRUE); 
cmd_gripe(Ieftover_argc,  leftoverargv); 
gcsjnode  =  CTDB_MODE_SIMNET ; 
reader_ihit(0); 

coordsettccawareness(TRUE); 
gcs_init(coord_get_datum_info); 
gcs_set_ceIl_awareness(GCS_MULTI_CELL); 
ctdb_cmplr_init(); 
num_correction_files  =  0; 

ctdbin  =  ctdb_cmplr_read(options.input_ctdb.value, 
ctdb_set_ti!ing_mode(FALSE); 
ctdbJntern_print_description(ctdb_in); 
ctdbjn->max_z  =  -999999.0; 
ctdb_in->min_z  =  999999.0; 
look_for_features(patch_water_state); 
for  (iy=0;  iy<=ctdb_in->max_y_post;  iy++) 

{ 


&input_format); 

j  This  loop  reads  the  elevations  from  the 
j  original  terrain  database ,  specified  by 
\  ctdbjn.  These  elevations  are  stored  in  the 
j  array  original  The  soil  for  each  grid  post  is 
\  initially  set  at  a  constant  value .  The  soil  is 
j  modified  by  transform _posts. 


This  section  processes  the  command  line 
argument  and  initializes  coordinate  system 
I  routines . 
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loc_y  =  iy*ctdb_in->incrD; 

for  (ix=0;  ix<=ctdb_in->max_xjpost;  ix++) 

{ 

Ioc_x  =  ix*ctdb_in->incrD; 

original[ix][iy]  =  ctdbJntern_get_ground_elevation 
(ctdb_in,loc_x,loc_y,  FAVORITES); 
if  ( original[ix][iy]  >  ctdb_in->max_z) 
ctdb_in->max_z  =  original[ix][iy]; 
if  (  original[ix][iy]  <  ctdb__in->min_z) 
ctdb_in->min_z  =  original  [ix]  [iy] ; 
my_soil[ix][iy]  =  131074; 

} 

} 


hgt_max  =  ctdb_in->max_z  -  ctdbjn->min_z; 
output_format  =  CTDBJFILEJFORMAT; 
dbase_type  =  CTDB  J3RIDDED; 
tcc  =  coord_define_tcc(COORD_UTM_NE, 

ctdb_in->origin_northing,ctdbJn->origin_easting, 
ctdb_in->origin_zone_number,ctdb_in->origin_zoneJetter, 
ctdb_in->datum,(ctdb_in->max_x  -  ctdb_in->min_x), 
(ctdb_in->max_y  -  ctdb_in->min_y)); 
getJIJ>ounds(COORDJTCC,  GCSJLLEGAL__CELL,  tcc, 
(float64)ctdbJn->min_x,  (float64)ctdbJn->min_y, 
(float64)ctdb_in->max_x,  (float64)ctdb_in->max_y, 
&indb_west_long,  &indb_south_lat, 

&indb_east_long,  &indb_north_lat); 
ncells  =  START_CELLS; 

if(!cells)  cells  =  (int32  *)STDALLOC(ncells  *  sizeof(int32)); 
if  (!gcs_extent_to_cell_list(indb_south  Jat,  indbwestlong, 
indb_northJat,  indb_east_Iong, 

&ncells,  cells)) 

{ 

cells  =  (int32  *)STDREALLOC(ceIls,  ncells  *  sizeof(int32)); 

} 

cells[0]  =  GCS_ILLEGAL_CELL; 
ncells  =  1 ; 
curcell  =  cellsfO]; 

dbname  =  ctdb_generate__tdb_filename(CTDB_MODEjSIMNET, 
options.tdb_path.value, 
options.output_base_name.  value, 
output  format,  GCS_ILLEGAL_CELL,  FALSE); 
header_assigned  =  TRUE; 
feature_offset  =  0; 
bzero((void*)&hdr,  sizeof(hdr)); 
my_hdr.max_x_post  =  2.0*ctdb_in->max_x_post; 
my_hdr.max_y_post  =  2.0*ctdb_in->max_y_post; 
my_hdr.pages_wide  =  my_hdr.max_x_post/32  +  1; 
my_hdr.pages_high  =  my_hdr.max_y_post/32  +  1; 
hdr.pages_wide  =  my_hdr.pages_wide; 
hdr.pages_high  =  my_hdr.pages_high; 
ctdb_in  ->  pages_wide  =  my_hdr.pages_wide; 
ctdb Jn  ->  pages_high  =  my_hdr.pages_high; 
ctdb_intern_pack_header(ctdb_in,  &hdr); 
hdr.format  =  output_format; 
time((time_t  *)&tloc); 


Myjidr  stores  the  information  needed  to 
increase  the  grid  size  for  the  output 
database ,  ctdb_out.  In  this  particular 
case ,  the  number  of  grid  posts  is  doubled . 
The  page  variables  specify  how  many 
pages  (a  32  x  32  array  of  integers )  are 
needed  to  store  the  final  database .  A 
patch  is  a  4x4  array  of  posts  used  to 
streamline  feature  lookup  routines . 


This  program  uses  only  one  cell 


This  section  initializes  coordinate 
transformation  routines . 
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strcpy(hdr.date,  ctime((time_t  *)&tloc)); 
hdr.num_linear_modeIs  =  CTDBJMAX_CANOPY_MODELS; 
sprintf(hdr.name,options.outputJ)ase_name.value); 
ctdbJntern_unpack_header(&ctdb_out,  &hdr); 
ctdb_out.eIevations  =  elev_ptr; 
ctdb_out.soils  =  soil_ptr; 

my_hdr.min_x  =  0.0;  my_hdr.min_y  =  0.0;  myjidr.min_z  =  0 
my_hdr.max_x  =  ctdb__in->max_x; 
my_hdr.max_y  =  ctdb__in->max_y; 
my_hdr.max_z  =  ctdb_in->max_z; 
myjidr.patches_wide  =  (my_hdr.pages_wide*32  +  2)/4; 
my_hdr.patches_high  =  (myjidr.pages_high*32  +  2)/4; 
my_hdr.incrD  = 

(my_hdr.max_x  -  my_hdr.min_x)/ 

(my_hdr.max_x__post  +  0.0); 
myjidr.  inv_incr  =  1.0/myJidr.incrD; 
myjidr.patchjncr  =  4; 

my_hdr.patch_inv_incr  =  1.0/(float)  (myjidr.patchjncr); 
myjidr.patch JncrD  =  my_hdr.patch_incr; 
my_hdr.incr  =  (float)  (myJidr.incrD  +  0.5); 
min_patch_x  =  0;  min_patch_y  =  0; 
max_patch_x  =  my_hdr.patches_wide; 
max_patch_y  =  myjhdr.patchesjiigh; 
ctdb_out.min_x  =  myjidr. min_x; 
ctdb_out.min_y  =  myjidr.min_y; 
ctdb_out.min_z  =  myjidr.min_z; 
ctdb_out.max_x  =  myjidr.max_x; 
ctdb_out.max_y  =  myjidr.max_y; 
ctdb_out.max_z  =  myjidr.maxjz; 
ctdb_out.pages_wide  =  my_hdr.pages_wide; 
ctdb_out.pagesJiigh  =  my_hdr.pages_high; 
ctdb_out.max_x_post  “  myJidr.max_x_post; 
ctdb_out.max_y_post  =  my_hdr.max_y_post; 
ctdb_out.patches_wide  =  myjidr.patches__wide; 
ctdb_out.patches_high  =  myjhdr.patchesjiigh; 
ctdb_out.incr  =  my_hdr.incr; 
ctdb_out.invJncr  =  myjidr.  in  vjncr; 
ctdb_out.incrD  =  myJidr.incrD; 
ctdb_out.patch_incr  =  myjidr.patchjncr; 
ctdb_out.patch  JnvJncr  -  my  Jidr.patchJnvJncr; 
ctdb_out. patch  JncrD  =  myjidr. patch  JncrD; 
ctdb_out.soilJables  =  (int32  **) 

ctdb_alloc(NULL,  ctdb_in->num_soil_tables  * 
sizeof(int32  *), "Compiler  rep.  soil  table”); 
ctdb_out . so i 1  Jab le_storage  =  (int8  *) 

ctdb_a!loc(NULL,  ctdb_in->soil_table__size, "Compiler  rep.  soil  table  data"); 
bcopy((void*)ctdbjn->soiljable_storage,(void*)ctdbj)utsoiljable_storage, 
ctdbJn->soiljable_size); 

ctdb_out.soil_tables[0]  =(int32  *)ctdbj>ut.soiljable_storage; 
for(i=l;i<ctdbJn->num_soilJables;i++) 

{ 

ctdb_out.soil_tables[i]  =  (int32  *  )(ctdb jwt .  so  i l_tab  1  e_storage  + 

((int8  *)ctdbJn->soilJables[i]  -  ctdbjn->soiljable_storage)); 

} 

ctdb_out.pat_table_ptr  =  (CTDB_PAT_TABLE_CMP  *) 
ctdb_a!loc(NULL,  sizeof(CTDB_PAT_TABLE_CMP), 


The  soil  table  is  copied  from  the  old 
terrain  database  to  the  new  database. 


These  variables  in  the  my_hdr 
structure  are  recalculated  based  on 
the  array  size  specified  by 
myjidr.postjiigh  and 
my Jidr. posts jwide.  The  structure  for 
the  output  database ,  ctdbjout ,  is 
updated  with  the  new  information 
stored  in  the  structure  myjidr. 


.0; 


Note:  myjidr.incr  is  an  integer 
whereas  myJidr.incrD  is  afloat64. 
It  is  important  that  myjidr.  incrD  is 
calculated  from  myjidr.  incr.  Both 
variables  are  used  by  the  ModSAF 
terrain  routines. 
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"Compiler  Polygon  Attribute  Table"); 
ctdb_out.pat_table_ptr->num_columns  =  ctdb_in->pat_tabIe_ptr->num_columns ; 
ctdb_out.pat_table_ptr**>alloced_columns  =  ctdb_in->pat_tabIe_j)tr->num_columns; 
ctdb_out.patJablejtr->num_entries  =  ctdb_in->pat_table_ptr->num__entries; 
ctdb_out.pat_table_ptr->alloced_entries  =  ctdb_in->pat_table__ptr->alloced_entries; 
ctdb_out.pat_table_ptr->columns  =  (CTDB_PAT_COLUMN_HDR_CMP  *) 
ctdb_alloc(NULL,  ctdb_put.patJable_ptr->aIloced_columns  * 
sizeof(CTDB_PAT_COLUMNJIDR_CMP), 

"Compiler  PAT  column  headers"); 
for(i“0;  i<ctdb_out.pat_table_ptr->num_columns;  i++) 

{ 

ctdb_out.pat_table_ptr->columnsti].facc  = 
ctdbJn->pat_table_ptr->columns[i].facc; 
ctdb_out.pat_table__ptr->columns  [i]  .data  - 
(CTDB_PAT_DATA_CMP  *)ctdb_alloc(NULL, 
ctdb_out.patJable_ptr->alloced_entries  * 
sizeof(CTDB_PATJ)ATA_CMP),  "Compiler  PAT  data"); 
bcopy((void*)ctdb_in->pat_table_ptr->columns[i].data, 
(void*)ctdb_out.pat_tablej)tr->columns[i].data, 
ctdb_in->pat_table_ptr->num_entries  * 
sizeof(CTDB_PAT_DATA_CMP)); 

} 

ctdb_out.water_chars  =  (CTDB  WATER  CHARS  CMP  *) 
ctdb_alloc(NULL,  ctdb_in->num_water_chars  * 
sizeof(CTDB__WATER_CHARS_CMP),"Cmp.  rep.  water  chars"); 
bcopy((void*)ctdbJn->water_chars,  (void*)ctdb_out.water_chars? 

ctdbjn->num_water_chars  *sizeof(CTDB_WATER_CHARS_CMP)) ; 
ctdb_outwater__chars[0].tidal_zonejndex  =  1; 
ctdb_out.tidal_zones  =  (CTDB_TIDAL_ZONE_CMP  *) 
ctdb_alIoc(NULL,  ctdb_in->num_tidal_zones  * 
sizeof(CTDB_TIDAL_ZONE_CMP),  ’ 

"Cmp.  rep.  water  chars"); 

bcopy((void*)ctdbJn->tidal_zones,  (vo  id  *  )ctdb_out.tidal_zones, 
ctdbJn->numjidal_zones  *sizeof(CTDBJTIDAL_ZONE_CMP)); 
ctdb_out.volume_models  =  NULL; 
ctdb_out.mes_voIume_models  =  NULL; 
ctdb_out.origin_x=-(ctdb_out.min_x+ctdb_out.max_x)/2; 
ctdb_out.origin__y=-(ctdb_out.min_y+ctdb_out.max_y)/2; 
ctdb_out.cell_id  =  GCS_ILLEGAL_CELL; 
ctdb__out.gcs_mode  =  gcs_mode; 
ctdb_out.west_long  =  indb_west Jong; 
ctdb_out.southJat  =  indb_south_lat; 
ctdb_out.east Jong  =  indb_eastJong; 
ctdb_out.north Jat  =  indb_north  Jat; 
ctdb_out.eIevations  =transform_posts(ctdbJn->elevations); 
ctdbjn->pages_wide  =  saved_pages_wide; 
ctdbjn->pagesjiigh  =  saved_pagesjiigh; 
patch_water_state  =  (int32  *) 

ctdb_alIoc(NULL,  (ctdb_out.patches_wide)*  (ctdb_out.patchesJiigh)* 
sizeof(int32),  "Patch  water  state  array"); 
bzero((void*)patch_water_state,  (ctdb_out.patches_wide)  * 

( ctdb_out. patches Jiigh)*  sizeof(int32)); 
encode_abstract_features(patch_water_state); 
add_subcanopies  (patch_water_state); 
hdr.num_features  =  encode_physicaI_features(output_format, 
dbname,  dbasejype,&hdr.numJinear_models, 


Encode jxbstract  features  adds  canopies 
and  other  polygonal  features  to  the 
database .  Encode jphysical features 
adds  roads ,  rivers,  trees  and  buildings 
to  the  database . 


Transform joosts  writes  the  new 
elevation  posts  to  ctdbjout 


The  water  characteristics  and  the 
tidal  zone  information  are  passed 
directly  from  ctdbjn  to  ctdbjout. 


The  polygonal  pattern  table  is 
transferred  from  ctdbjn  to  ctdbjout. 
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&hdr.num_aggregate_models,  patch_water_state) ; 
ctdb_encode_ab$tract(&ctdb_out,  &quad_root); 
hdr.num_quad_data  =  ctdb_outnum_quad_data; 
collect_nodes_and_edges(&hdr); 

ctdb_cmplr_write(&ctdb_out,  dbname,  outputjormat); 
ctdb_j>rint_feature_stats(); 

if  (ctdb_out.features) 

munmap(ctdb_out.features, 

hdr.num_features  *  sizeof(CTDB_FEATURE_DATA_CMP)); 
if  (patch_water_state) 

STDDEALLOC(patch_water_state); 

} 

return  0; 


A.3  Reorder_patch_features 

static  void  reorder_patch_features  (offsets,  sizes,  dbname,  num_features) 

int32  *offsets; 
int32  *sizes; 
char  *dbname; 
int32  num_features; 

{ 

FILE  *fp; 

char  temp_fhame[255]; 

int32  patch_x,  patch_y,  patch_number,  offset  =  0; 
sprintf(temp_fhame,  "%s.features",dbname); 
l\vrite(ctdb_out. features,  sizeof(CTDB_FEATURE_DATA_CMP), 
num_features,  fp); 
fclose(fp); 

for(patch_y  =  0;  patch_y  <  ctdb_out.patches_high;  patch_y++) 

{ 

for(patch_x  =  0;  patch_x  <  ctdb_out.patches_wide;  patch_x++) 

{ 

patch_number  =  (patch_y  *  ctdb_out.patches_wide)  +  patch_x; 
if(offsets[patch_number]  !=  NO_FEATURES) 

{ 

fseek(fp,  offsets[patch_number]  * 
sizeof(CTDB_FEATURE_DATA_CMP),0); 
fread(ctdb_out.features  +  offset, 
sizeof(CTDB_FEATURE_DATA_CMP), 
sizes[patch_number],  fp); 

ctdb_set_patch_feature_start(&ctdb_out,  patchnumber,  offset); 
offset  +=  sizes[patch_number]; 

} 

} 

} 

fclose(fp); 

unlink(temp_fhame); 

} 

A.4  Encode_physical_features 

static  int32  encode_physical_features  (format,  dbname,  dbase_type, 
num_linear_mode!s,  num_aggregate_models,  patch_water_state) 
int32  format; 


Reoder jpatchjeatures  is 
unchanged  from  the  original 
program  recompile,  c  in  the 
ModSAF  3.0  source  code. 
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char  *dbname; 

CTDBDBASETYPE  dbase_type; 
int32  *num_linear_models; 
int32  *num_aggregate_models; 
int32  *patch_water_state; 

{ 

int32  patches; 

int32  patch_meters,  oldj5atch_meters,new_patch_meters; 
int32  patch_x,  patch_y,  new_patch_x,  new_patch_y; 
int32  old_patch_x,  old_patch_y; 
int32  added  =  0; 

CTDB_LAID_LINEAR  ^linear,  *all_linear; 

char  temp_fiiame[CTDB_FILENAME_DIMENSION_CMP]; 

int32  ij,gcs_patch_num; 

static  CTDB_LAID_LINEAR  local; 

CMP_REF  tz_refs[CMP_MAX_REF] ; 
int8  *post_array  =  NULL; 
int32  *offsets,  *sizes; 

trans_patches  =  (TRANSITIONJPATCHINFO  *); 

ctdb_al»oc(NULL,sizeof(TRANSITION_PATCH_INFO)* 
ctdb_out.patches_wide  *  ctdb_out.patches_high,  "trans  patches")? 
bzero((void*)trans_patches,  sizeof(TRANSITION_PATCH_INFO)  * 
ctdb_out.patches_wide  *  ctdb_out.patches_high); 
offsets  =  (int32  *)ctdb_alIoc(NULL,  sizeof(int32)  * 

ctdb_out.patches_wide  *ctdb_out.patches_high,  "offsets"); 
for(i  =  0;  i  <  ctdb_out.patches_high  *  ctdb_out.patches_wide;  i++) 
offsets[i]  =  NO_FEATURES; 
sizes  =  (int32  *)ctdb_alloc(NULL,  sizeof(int32)  * 

ctdb_out.patches_wide  *ctdb_out.patches_high,  "sizes"); 
bzero((void*)sizes,  sizeof(int32)  *  ctdb_outpatches_wide  * 
ctdbout.patcheshigh); 

ctdb_zero_out_patch_bufferO? 

sprintf(temp_ffiame,  "%s.te",dbname); 

patches  =  ctdb_out.patches_wide  *  ctdb_out.patches_high; 

if(!ctdb_out.patch_groups) 

{ 

ctdb_out.patch_groups  =(CTDB_PATCH_GROUP_CMP  *)STDALLOC( 
sizeof(CTDB_PATCH_GROUP_CMP)* 
(patches/CTDB_PATCHES_PER_GROUP_CMP)); 
bzero((void*)ctdb_out.patch_groups,sizeof(CTDB_PATCH_GROUP_CMP)* 
(patches/CTDB_PATCHES_PER_GROUP_CMP)); 

} 

if(!ctdb_out.linear_models) 

{ 

ctdb_out.linear_models  =  (CTDB_LINEAR_MODEL_CMP  *) 
STDALLOC(MAX_NUM_LINEAR_MODELS  * 
sizeof(CTDB_LINEAR_MODEL_CMP)); 
bzero((void*)ctdb_out.linear_models,  MAX_NUM_LINEAR_MODELS  * 
sizeof(CTDB_LINEAR_MODEL_CMP)); 

} 

if( !  ctdb_out.aggregate_models) 

{ 

ctdb  out.aggregate  models  =  (CTDB_AGGREGATE_MODEL_CMP  **) 
STDALLOC(MAX_NUM_AGGREGATE_MODELS  * 
sizeof(CTDB_AGGREGATE_MODEL_CMP  *)); 
bzero((void*)ctdb_out.aggregate_models, 
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MAX_NUM_AGGREGATE_MODELS  * 
sizeof(CTDB_AGGREGATE_MODEL_CMP  *)); 

} 

old_patch_meters  =  ctdb_in->incr  *  ctdb_in->patch_incr; 
new_patch_meters  =  ctdb_out.incr  *  ctdb_out.patch_incr; 

for  (new_patch_y  =  0;  new_patch_y  <  ctdb_out.patches_high;  new_patch _y++,  putchar('.'),fflush(stdout)) 

{ 

for  (new_patch_x  =  0;  new_patch_x  <  ctdb_out.patches_wide;  new_patch_x++) 

{ 

old_patch_x  =  xx/old__patch_meters; 
old_patch_y  =  yy/oldjpatch_meters; 
patch_x  =  oId_patch_x; 
patch_y  =  old_patch_y; 

gcs_patch_num  =  patch_y  *  ctdb_in->patches_wide  +  patch_x; 
ctdb_initjpatch_bufferO; 

collect_canopies(new_patch_x,new_patch_y,patch_x,  patch  _y); 
collect_building$  (new_patch_x,new_patch_y,  old_patch_x,oId_patch_y); 
collect_trees  (new_patch_x,new_patch_y,  old_patch_x,old_patch_y); 
collectjinears  (new_patch_x,new_patch_y,  old_patch_x,old_patch  _y); 
if  (( old_patch_y  <  ctdb_in->patches_high)  &&  ( old_patch_x  >  0)) 

{ 

coiiect_building$  (new_patch_x,new_patch_y,  old_patch_x-l,old_patch_y+l); 
co!lect_trees  (new_patch_x,newjpatch_y,  oldjpatch_x-l,old_patch_y+l); 
coIlectJinears(new_patch_x,new_patch_y,  old_patch_x-l,oId_patch_y+l); 

} 

if  (  old_patch_y  <  ctdb_in->patches_high) 

{ 

collect_buildings  (new_patch_x,new_patch_y,  old_patch_x,old_patch_y+l); 
collect_trees  (new_patch_x,new_patch_y,  old_patch_x,old_patch _y+l); 
collectjinears  (new_patch_x,newjpatch_y,  old_patch_x,old_patch_y+l); 

} 

if  (( oId_patch_y  <  ctdb_in->patches_high)  &&  (  old_patch_x  <  ctdb_in->patches_wide)) 

{ 

collectjbuildings  (new_patch_x,new_patch_y,  oId_patch_x+l,old_patch_y+l); 
col!ect_trees  (new_patch_x,new_patch_y,  old_patch_x+l,old_patch_y+l) 
co!lect_linear  (new_patch_x,new_patch_y,  old_patch_x+l,old_patch_y+l); 

} 

if  (  old_patch_x  <  ctdb_in->patches_wide) 

{ 

collectjbuildings  (new_patch_x,new_patch_y,  old_patch_x+l,old_patch_y); 
colIect_trees(new_patch_x,new_patch_y,  oIdjpatch_x+l,old_patch  _y); 
collectjinear  (new_patch_x,new_patch_y,  old_patch_x+l,old_patch _y); 

} 

if  (( old_patch_y  >  0)  &&  ( old_patch_x  <  ctdbjn->patches_wide)) 

{ 

collectjbuildings  (new_patch_x,new_patch_y,  old_patch_x+l,old_patch_y-l); 
collecttrees  (new_patch_x,new_patch_y,  old_patch_x+l,old _patch_y-l); 
collectjinears  (new_patch_x,new_patch_y,old_patch_x+l,old_patch_y-l); 

} 

if  ( old__patch_y  >  0) 

{ 

collectjbuildings  (new_patch_x,new_patch_y,  old_patch_x,old_patch_y-l); 
collect Jrees  (new_patch_x,new_patch_y,  old_patch_x,old_patch_y-l); 
collectjinears  (newjpatch_x,new_patch_y,  old_patch_x,old_patch_y-l); 

} 

if  ((  old_patch_x  >  0)  &&  (  old_patch_y  >  0) 


The  new  patches  may  overlap  more  than  one  of  the  original 
patches,  so  features  must  be  collected  for  each  of  the  nine 
original  patches  that  might  intersect  the  new  patch. 
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collect_buildings  (new_patch_x,new_patch_y,  old_patch_x-l,old_patch_y  - 1); 
collect_trees  (new_patch_x,new_patch_y,  old_patch_x-l,old_patch_y  - 1); 
colIect_linears(new_patch_x,new_patch_y,  old_patch_x-l,old_patch_y  - 1); 

} 

if  (  old_patch_x  >  0) 

{ 

collect_buildings(new_patch_x,new_patch_y,  old_patch_x-l,old_patch_y); 
collect_trees  (new_patch_x,new_patch_y,  old_patch_x-l,oId_patch_y); 
collect_Iinears(new_patch_x,new_patch_y,  old_patch_x-l,old_patch_y); 

} 

added  =  ctdb_store_patch_buffer(&ctdb_out, 

new_patch_x  +  offset_x,  new_patch_y  +  offset_y, 
ctdb_out.features+feature_offset,  ctdb_out.linear_models, 
&ctdb_out.num_linear_modeIs,  ctdb_out.aggregate_models, 

&ctdb_out.num_aggregate_models,complete_micro_poly,  _ 

&num_te,  transjpatches,  &trans_end, 

&max_elev,  format,  dbase_type,  FALSE, 

FALSE,  fptr,  FALSE); 
if(added) 

{ 

offsets[(new_patch_y  +  offset_y)  *  ctdb_out.patches_wide  + 
new_patch_x  +  offset_x]  =  feature_offset; 
sizes [(new_patch_y  +  offset_y)  *  ctdb_out.patches_wide  + 
new_patch_x  +  offset_x]  =  added; 

} 

feature_offset  +=  added; 

} 

} 

if(ctdb_in->features)  STDDEALLOC(ctdb_in->features); 
if  (ctdb_in->elevations)  STDDEALLOC(ctdb_in->elevations); 
if  (ctdb_in->soils)  STDDEALLOC(ctdb_in->soils); 
if  (post_array)  STDDEALLOC(post_array); 
reorder_patch_features(offsets,  sizes,  dbname,  feature_offset); 

STDDEALLOC(offsets); 

STDDEALLOC(sizes); 
fclose(fptr); 

ctdb_out.n_te_bytes  =  total_tes  *  sizeof(TERRAIN_ELEMENT_16); 

*num_linear_models  =  ctdb_out.num_linear_models; 

*num_aggregate_models  =  ctdb_out.num_aggregate_models; 
return  feature_offset; 

} 

A.5  Look_for_features 
static  void  Iook_for_features  (patch_water_state) 
int32  *patch_water_state; 

{ 

CTDB_SEARCH_SPACE_PTR_CMP  searchspace; 
int32  i,  j,ix,iy,ii,max_steps,flag; 
int32  count, code,  data_size; 
float32  * verts,  *1  verts  =  NULL; 

float32  xl  ,y  1  ,x2,y2,cx,cy,delx,dely,slope,step,my_sign; 

CTDBABSTRACTDATACMP  *data; 

I  int32  num_lverts  =  0,  *soil_table,soil; 
float64  meters_per_patch; 

meters_per_patch  =  ctdb_in->incrD  *  ctdb_in->patch_incrD; 


Ctdb_store _patch_buffer  writes  all 
the  physical  features  to  ctdbjout. 
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step  =  ctdb_in->incrD/10.0; 
search_space  = 

ctdb_intern_create_search_space(ctdb_in,  ctdb_in->min_x,  ctdb_in->min_y, 
ctdb_in->max_x,  ctdbjn->max_y,  CTDBJBY_PATCH); 
while  (ctdb_intem_next_quad_patch(search_space)) 

{ 

for  (code=l  ;code<=24;code-H-) 
while  (count  =  ctdbjntern_next_abstract(search_space, 

{ 

num_lverts  =  count; 

if  (code  =  CTDB_ABSTRACT_CANOPY_CMP) 

{ 

for  (i=2 ;  i<(2  *  count) ;  i+=2) 

{ 

xl  =  verts[i-2];  yl  =  verts[i-l]; 
x2  =  verts[i];  y2  =  verts[i+l]; 
delx  =  x2-xl ;  dely  =  y2-y  1 ; 
if  ((delx*delx)  >  (dely*dely)) 

{ 

slope  =  dely/delx; 
my_sign=  1.0; 
if  (delx  <  0)  my_sign  =  -1.0; 
max_steps  =  sqrt(delx*delx)/(2.0*mat_step); 
for  (ii  =  0;  ii  <  max_steps;ii++) 

{ 

cx  =  xl  +  ii*step*my_sign; 
cy  =  slope*(cx  -  xl)  +  y  1 ; 
ix  =  (float)  (cx/mat__step); 
y  =  (float)  (cy/mat_step); 
canopy  [ix][iy]  =  1; 

} 

} 

else 

{ 

slope  =  delx/dely; 
my_sign  =  1.0; 
if  (dely  <  0)  my_sign  =  -1.0; 
max_steps  =  sqrt(dely*dely)/(2.0*mat_step); 
for  (ii  =  0;  ii  <  max_steps;ii++) 

{ 

cy  =  yl  +  ii*step*my_sign; 
cx  =  slope*(cy  -  yl)  +  xl ; 
ix  =  (float)  (cx/mat_step); 
iy  =  (float)  (cy/mat_step); 
canopy[ix][iy]  =  1; 

} 

} 

} 

} 

} 

} 

for  (iy  =  0;iy  <  5500;iy++) 

{ 

flag  =  0; 

if  ( canopy[0][iy]  =  1)  flag  =  1; 
for  (ix  =  1;  ix  <  5500;ix++) 


code,  Averts,  &data,  &data_size)) 

Look Jor features  finds  all  the 
canopies  and  lakes  on  the  original 
database  so  that  the  soil  type  for 
the  underlying  posts  can  be 
adjusted.  The  array  canopy  is  a 
high  resolution  map  of  the  location 
of  the  canopies. 
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{ 

if  ((canopy[ix][iy]  —  0)  &&  (flag  ==  1)) 
canopy[ix][iy]  =2; 

if  ((canopy  [ix][iy]  =  1)  &&  (flag  ==  0)) 
flag=  1; 

else  if  ((canopy[ix][iy]  =  1)  &&  (flag  =  1)) 
flag  =  0; 

} 

} 

ctdb_intern_destroy_search_space(search_space); 
if  (num_lverts  >  0)  free(lverts); 

} 

A.6  Encode_abstract_features 
static  void  encode_abstract_features  (patch_water_state) 
int32  *patch_water_state; 

{ 

CTDB_SEARCH_SPACE_PTR_CMP  search_space; 

int32  i,  j; 

int32  count; 

int32  code; 

float32  *verts; 

CTDB_ABSTRACT_DATA_CMP  *data; 

int32  data_size; 

float32  *1  verts  =  NULL; 

int32  num_lverts  =  0; 

int32  *soil_table; 

int32  soil; 

float64  meters_per_patch; 

meters_per_patch  =  ctdb_out.incrD  *  ctdb_out.patch_incrD; 
searchspace  =  ctdb_intern_create_search_space(ctdb_in, 
ctdb_outmin_x,  ctdb_outmin_y, 
ctdb_out.max_x,  ctdb_out.max _y, 
ctdb_BY_PATCH); 

while  (ctdb_intern_next_quad_patch(search_space)) 

{ 

for  (code=l  ;code<CTDB_ABSTRACT_MAX;code++) 

while  (count  =  ctdb_intern_next_abstract(search_space,  code,&verts,  &data,  &data_size)) 

{ 

if  (count  >  num_lverts) 

Averts  =  (float32*)STDKEALLOC(lverts,2*count*sizeof(float32)); 
num_lverts  =  count; 

} 

for(i=0;i<(2*count);i+=2) 

{ 

lverts[i]  =  verts[i]  +  offset_x_meters; 
lverts[i+l]  =  verts[i+l]  +  offset__y_meters; 

} 

if  (code  ==  CTDB  ABSTRACT  CANOPY  CMP) 
data->canopy.  impenetrable  =  1; 

ctdb_store_abstract(&ctdb_out,  &quad_root,code,  count,  Iverts, 
CTDB_ABSTRACTJDATA_SIZE(data_size), 

(int32  *)data,  min_patch_x  *  meters_per_patch, 
min_patch_y  *  meters_per_patch,max_patch_x  *  meters_per_patch, 
max_patch_y  *  meters_per_patch); 


Encode _abstract  tranfers  the  data  for  lakes, 
canopies,  and  other  abstract  features  from 
ctdbjn  to  ctdb_out.  Note  that  the  polygons 
do  not  depend  on  patch  size  so  they  do  not 
need  to  be  adjusted  for  the  new  database. 
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} 


} 

ctdb_intern_destroy_search_space(search_space); 
if  (num_lverts  >  0)  free(lverts); 

} 

A.  7  CoIlect_nodes_and_edges 
static  void  co!lect_nodes_and_edges(hdr) 
CTDB_FILE_HEADER  *hdr; 

{ 

ctdb_create_networks(&ctdb_out,  output_format); 
hdr->num_nodes  =  ctdb_out.num_nodes; 
hdr->node_size  =  ctdb_out.node_size; 
hdr->num_edges  =  ctdb_out.num_edges; 
hdr->edge_size  =  ctdb_out.edge_size; 

} 


A.8  Transform  jposts 

static  int32  *transform_posts() 

{ 

float64  x,y,dx,dy»zz,z00,zl  0,z01,zl  1 ; 
float64  xd,yd,loc_x,  loc__y; 

float64  page_meters,  post_meters,inv_post_meters; 
uint32  p,  s; 

int32  *new_posts,  *new_soils; 
int32  norm_x,  norm_y  .ix.iy; 
int32  post_index,page,  index; 

new_posts  =  (int32  *)ctdb_alloc(NULL,  ctdb_in->pages_wide  * 
ctdb_in->pages_high  *  PAGE_SIZE,"New  elevation  posts"); 
new_soils  =  (int32  *)ctdb_alloc(NULL,  ctdb_in->pages_wide  * 
ctdb_in->pages_high*  PAGE_SIZE,  "New  soil  indices"); 
post_meters  =  ctdb_out.incrD; 
page_meters  =  POSTS_PER_SIDE  *  post_meters; 
inv_post_meters  =  1.0/post_meters; 
add_hills_to_patches(); 
for  (y=0;  y<=ctdb_out.max_y;  y+=page_meters) 
for  (x=0;  x<=ctdb_out.max_x;  x+=page_meters) 
for  (dy=0.0;  dy<page_meters;  dy+=post_meters) 
for  (dx=0.0;  dx<page_meters;  dx+=post_meters) 

{ 

loc_x  =  x+dx; 
loc_y  =  y+dy; 

ix  =  (float)  (loc_x/ctdb_in->incrD); 
iy  =  (float)  (loc_y/ctdb_in->incrD); 
xd  =  (loc_x  -  (float)  (ix*ctdb_in->incrD)); 
yd  =  (loc_y  -  (float)  (iy*ctdb_in->incrD)); 
zOO  =  xd*(original[ix  +l][iy]  -  original[ix][iy])/ctdb_in->incrD+ 

yd*  (original  [ix][iy  +  1]  -  original[ix][iy])/ctdb_in->incrD  +  original  [ix][iy]; 
zl  1  =  xd*  (original[ix][iy+l]  -  original[ix+l][iy+l])/ctdb_in->incrD+  yd* 

(original[ix+l][iy]  -  original[ix+l][iy+l])/ctdb_in->incrD  +  original[ix+l][iy+l]; 
norm  x  =  (int32)POST_X(loc_x  *  inv _post_meters); 
norm_j/  =  (int32)POST_Y(loc_y  *  inv_post_meters); 
zz=  (zOO  +zll)/2.0; 
zz  =  zz  +  add_vrt(loc_x,_locy); 


Addjiillsjo _patches  adds  vrt  hills  to 
each  patch  so  that  the  original  elevation 
data  are  not  altered.  Transform  jposts 
produces  the  elevation  grid for  ctdbjout 
from  the  original  elevation  grid  stored  in  j 
original  and  the  new  vrt  hills.  This 
routine  also  adjusts  the  soil  type,  based  j 
on  the  array  my_soil. 


Collect _nodes_and_edges  is  copied  from 
the  original  program  recompile,  c  in  the 
ModSAF  3.0  source  code. 
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p  =  BUILD  JPOSTJELEVATION  (zz,  ctdbjn->inv_fixed _point_basis); 
p  &=  ~POST_FLAGS; 
if  (canopy  [ix][iy]  =  1)  p  |=  POSTJTREES; 
s  =  my_soil[temp_x][temp_y]; 

post_index  =  ctdb_intern_get_post_index(ctdb_in,  norinjt,  norm_y); 
new_posts[post__index]  =  p; 
new_soils[post_index]  =  s; 

} 

ctdb_out.soils  =  new_soiIs; 
return  new__posts; 

} 

A.9  Add_hilIs_to_patches 

add_hills_to_patchesO; 

{ 

int  ix,iy,ih; 

float64  w,base_hgt,xx,yy,xO,yO,xl,yl; 
float64  d,d00,d01,dll,dl0; 
w  =  ctdb_in->patch_width; 
base_hgt  =  5.0; 

for  (iy  =0;iy<  ctdb_in->patches_high;iy++) 

{ 

for  (ix=0;ix  <  ctdb_in->patches_wide;ix-H-) 

{ 

hills_per_patch[ix][iy]  =  5.0*drand48(); 
for(ih=0;ih  <  hills_per_patch[ix][iy];ih++) 

{ 

hill[ihill].cx  =  xx=  w*drand48(); 
hill[ihill].cy  =  yy=  w*drand48(); 
hill[ihill].hgt  =  base_hgt*drand48(); 
hill[ihill]  .angle  =  1 80.0*drand48(); 
hill[ihill].wgty  =  0.75*drand48()  +  0.25; 
hill[ihill].power  =  5.0*drand48()  +  0.5; 
xO  =  (int)  (hill[ihill].cx/ctdb_in->incrD); 
yO  =  (int)  (hill[ihill].cy/ctdb_in->incrD); 
xl  =  xO  +  ctdb_in->incrD; 
yl  =  yO  +  ctdb_in->incrD; 
dOO  =  sqrt((xx-x0)  *  (xx-xO)  +(yy-yO)*(yy-yO)); 
dlO  =  sqrt((xx-xl)*(xx-xl)  +(yy-yO)*(yy-yO)); 
dl  l  =  sqrt((xx-xl)*(xx-xl)  +(yy-yl)*(yy-yl)); 
d01=  sqrt((xx-xO)  *  (xx-xO)  +(yy-yl)*(yy-yl)); 
d  =  dOO; 

if  (d  >  dlO)  d  =  dlO; 

if(d>  dll)d  =  dll; 

if  (d  >  dOl)  d  =  dOl; 

hill[ihill].wgtx  =  d*drand48()  +  0.5; 

hill[ihill].cx  =  hill[ihill].cx  +  ctdb_in->patch_width; 

hill[ihill].cy  =  hill[ihill].cy  +  ctdb_in->patch_width; 

} 

} 

} 

} 

} 

A.10  Compiete_micro_poIy 

static  void  complete_micro_poly(poly) 


Complete _micro _poly  is  copied  from  the 
original  program  recompile,  c  in  the  ModSAF 
3.0  source  code. 
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CTDB_MC_POLY  GON  *poly; 

{ 

int32  i; 

int32  norm_x,  norm_y; 
float64  midx,  midy; 
midx  =  midy  =  0.0; 
for  (i=0;i<3;i++) 

{ 

norm_x  =  poly->verts[i][X]  *  ctdb_out.inv_incr; 
norm_y  =  poly->verts[i][Y]  *  ctdb_out.inv_incr; 
if  (gcsjnode  =  CTDB_MODE_SIMNET) 

{ 

poly->verts[i][Z]  =ctdb_intern_lookup_post_elevation(ctdb_in,  norm_x,  norm_y); 

} 

else 

{ 

poly->verts[i][Z]  =ctdb_intern_Iookup_post_elevation(&ctdb_out,  norm_x,  norm_y); 

} 

midx  +=  norm_x; 
midy  +=  norm_y; 

} 

midx  /=  3.0; 
midy  /=  3.0; 

poly->soil  =  ctdb_intern_lookup_soil(ctdb_in,  midx,  midy); 

} 

A.11  Collect_building 

static  void  collect_buiIding  (new_patch_x,new_patch_y,old_patch_x,old_patch_y) 
int32  new_patch_x,  new_patch_y ,  old_patch_x,old_patch_y; 

{ 

CTDBVOLUMEPOLYGON  building; 
int32  n_verts,  i,  j,  offset; 

CTDB_FEATURE_DATA_CMP  *  feature,  *last; 

int32  patch_number  =  old_patch_x  +  old_patch_y  *  ctdb_in->patches_wide; 
int32  next,  fclass,x,y; 

int32  meters_per_patch  =  ctdb_in->incr  *  ctdb_in->patch_incr; 
int32  patch_min_x  =  old_patch_x  *  meters_perjpatch; 
int32  patch_min_y  =  old__patch_y  *  meters__per_patch; 
int32  patch_max_x  =  (old_patch_x+l)  *  meters_per_patch; 
int32  patch_max_y  =  (old_patch_y+l)  *  meters_per_patch; 
int32  intersection; 
offset  =  3; 

feature  =  ctdb_intem_lookup_patch(ctdb_in,  patch_number); 
if  (feature) 

{ 

last  =  feature  +  feature->hdr.size; 
feature++; 

for(;  feature  <  last;  feature  +=  next) 

{ 

n_verts  =  feature->info.vertices; 
fclass  =  feature->info.feature_class; 
if  (fclass  =  F  C_MICRO_CMP) 

{next  =  NEXT_MICRO_CMP(n_verts);continue; } 
else  if  (fclass  ==  FC_LINEAR_CMP) 

{next  =  NEXT_LINEAR_CMP(n_verts);continue; } 
else  if  (fclass  ==  FC_CANOPY_CMP) 


This  section  skips  all 
features  except  the 
buildings 
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} 


} 


{next  =  NEXT_CANOPY_CMP(n_verts);continue; } 
else  if  (fclass  =  FC_AGGREGATE_CMP) 

{next  =  NEXT_AGGREGATE_CMP(n_verts);continue; } 
else  if  (fclass  ==  FC_LAID_LINEAR_CMP) 

{next  =  NEXT_LAID_LINEAR(n_verts);  continue;  } 
else  if  (fclass  ==  FC_VOLUME_CMP) 
next  =  NEXT_VOLUME_CMP(n_verts); 
if  (n_verts) 

{ 

recompile_buildings++; 
building.n_verts  =  n_verts; 
for  (i=0;i<n_verts;i++) 

{ 


building. verts[i][X]  =  ctdb_in->incrD  * 

X  YO_TO_POST(feature  [offset+2  *i  ]  .xy.x,old_patch_x, 
ctdb_in->patch_incrD)  +  offset_x_meters; 
building.  verts[i][Y]  =  ctdb_in->incrD  * 

XYO_TO_POST(feature[offset+2*i].xy.y,old_patch_y, 
ctdb_in->patch_incrD)  +  offset_y_meters; 
building.verts[i][Z]  =  feature[offset+l+2*i].z; 

} 

if  ( n_verts  <  3) 

{ 

building.n_verts  =  3; 

building. verts[2][X]  =  10.0  +  ctdb_in->incrD  * 

XYO_TO_POST(feature[offset].xy.x,old _patch_x, 
ctdb_in->patch_incrD)  +  offsetxmeters; 
building.verts[2][Y]  =  ctdb_in->incrD  * 

XYO_TO_POST(feature[offset].xy.y,old_patch_y, 
ctdb_in->patch_incrD)  +  offset_y_meters; 
building.  verts[2][Z]  =  feature[offset+l].z; 

} 

building.type  =  feature[l].volume.type; 
building.reference  =  feature[l].volume.reference; 
ctdb_buffer_vol(&building); 


Volume  features  are  sets 
of  polygon  vertices  in 
three-dimensional  space. 
For  each  voulme  feature, 
building  stores  those 
vertices  that  intersect  the 
current  patch.  If  there  are 
fewer  than  three  vertices, 
an  extra  vertex  is  added  to 
complete  the  polygon. 


} 


A.12  CoIlect_trees 

static  void  collect_trees  (new_patch_x,new_patch_y,old_patch_x,old_patch_y) 
int32  newjpatch_x,  new_patch_y,old_patch,old__patch_y; 

{ 

CTDBJLINEAR  tree; 
int32  n_verts,  i,ix,iy; 

CTDB_FEATURE_DATA_CMP  ♦feature,  *last; 

int32  patch_number  =  old__patch_x  +  old_patch_y  *  ctdb_in->patches_wide; 
int32  next,  fclass; 

int32  my_max_x,  my_max_y,  my_min_x,  my_min_y; 
float64  xx,yy,meters_per_patch; 

my_min_x  =  new_patch_x*ctdb_out.patch_incr*ctdb_out.incr; 
my_min_y  =  new_patch_y*ctdb_out.patch_incr*ctdb_out.incr; 
my_max_x  =  (new_patch_x  +  l)*ctdb_out.patch_incr*ctdb_out.incr; 
my_max_y  =  (new_patch_y  +  l)*ctdb_out.patch_incr*ctdb_out.incr; 
meters__per_patch  =  ctdb_out.incrD  *  ctdb_out.patch_incrD; 
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feature  =  ctdb_intem_lookup_patch(ctdb_in,  patch_number); 
if  (feature) 

{ 

last  =  feature  +  feature->hdr.size; 

feature++; 

for(;  feature  <  last;  feature  +=  next) 

{ 

n_verts  =  feature->info.vertices; 
fclass  =  feature->info.feature_class; 
if  (fclass  ==  F  C_MICRO_CMP) 

{next  =  NEXT_MICRO_CMP(n_verts);continue; } 
else  if  (fclass  =  FC_VOLUME_CMP) 

{next  =  NEXT_VOLUME_CMP(n_verts);continue; } 
else  if  (fclass  ==  FC_CANOPY_CMP) 

{next  =  NEXT_CANOPY_CMP(n_verts);continue; } 
else  if  (fclass  =  FC_AGGREGATE_CMP) 

{next  =  NEXT_AGGREGATE_CMP(n_verts);continue;} 
else  if  (fclass  =  F C_L  AID_LINEAR_CMP) 

{next  =  NEXT_LAID_LINEAR_CMP(n_verts);continue; } 
else  if  (fclass  =  FC_LINEAR_CMP) 

next  =  NEXT_LINEAR_CMP(n_verts); 
if  (n_verts) 

{  “ 

recompile_trees-H-; 
tree.n_verts  =  n_verts; 
tree.desc.type  =  CTDB_FT_TREELINE_CMP; 
tree.desc.linear_data.treeline_info.foliage_height  = 
ctdb_in->incrD*  X  Y  0_T  0_P0ST  (ctd  b_in->linea  r_m  odels  [ 
feature{l].linear.reference).Iinear_data.treeline_info. 
foliageheight,  0,  ctdb_in->patch_incrD); 
tree.desc.linear_data.treeline_info.trunk_radius  = 

ctdb_in->incrD*XYO_TO_POST(ctdb_in->linear_models[ 
feature[l].Iinear.reference].linear_data.treeIine_info. 
trunk_radius,  0,  ctdb_in->patch_incrD); 
tree.desc.width  =  ctdb_out.incrD  *  XYO_TO_POST( 

ctdb_in->linear_models[feature[l].Iinear.referenceJ.width, 
0,  ctdb_in->patch_incrD); 

tree.desc.linear_data.treeline_info.fullness  = 

ctdb_in->linear_models[feature[l].linear.reference]. 

linear_data.treeline_info.fullness/ 

(float32)(l  «  CTDB_FULLNES S_BITS_CMP  - 1); 
tree.desc.linear_data.treeline_info.total_height  =  0.0; 
for  (i=0;i<n_verts;i++) 

{ 

tree.verts[i][X]  =  ctdb_in->incrD  * 

XYO_TO_POST(feature[2+2*i].xy.x,old_patch_x, 
ctdb_in->patch_incrD)  +  offset_x_meters; 
tree.verts[i][Y]  =  ctdb_in->incrD  * 

XYO_TO_POST(feature[2+2*i].xy.y,old_patch_y, 
ctdb_in->patch_incrD)  +  offset_y_meters; 
tree.verts[i][Z]  =  feature[3+2*i].z; 

} 

ctdb_buffer_linear(&tree); 
if  (n_verts  =1) 


Individual  trees  are  stored  for  the 
first  and  last  tree  in  a  tree  line. 


Tree  lines  are  linear  segments  with 
information  about  foliage  height, 
foliage  density,  tree  spacing,  and  trunk 
radius. 


This  section  skips  all 
features  except  the  tree  lines 
and  individual  trees. 
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{ 

tree.n_verts  =  1 ; 

tree.verts[l][X]  =  tree.verts[0][X]  +20.0; 
tree.verts[l][Y]  =  tree.verts[0][Y]; 
ctdb_buffer_Iinear(&tree); 

} 

if  ((n_verts  >  1)  &&  (feature[l].integer  &  M_TREE_FIRST_CMP)) 

{ 

ctdb_buffer_linear(&tree); 
tree.n_verts  =  1 ; 
ctdb_buffer_linear(&tree); 

} 

if  ((n_verts  >  1)  &&  (feature[l].integer  &  M_TREE_LAST_CMP)) 

{ 

tree.n_verts  =  1 ; 

tree.verts[0][X]  =tree.verts[n_verts-l][X]  +offset_x_meters; 
tree.verts[0][Y]  =  tree.verts[n_verts-l][Y]  +offset_y_meters; 
tree.verts[0][Z]  =  tree.verts[n_verts-l][Z]; 
ctdb_buffer_linear(&tree); 

} 

} 

} 

} 

} 

A.13  Collect_linear 

static  void  collect_linear  (new_patch_x,new_patch_y,oldjpatch_x,old_patch_y) 
int32  new_patch_x,new_patch_y,old_patch_x,old_patch _y; 

{ 

CTDB_LAID_LINEAR  line; 
float64  verts[CTDB_VERTICES_MAX]  [2] ; 
float64  width, xx, yy,xl,yl,x2,y2; 
float64  xm,ym,xM,yM,del_x,  del_y; 
int32  soil,n_verts„k,j; 

CTDBFEATUREDATACMP  ^feature,  *last; 

int32  patch_number  =  old_patch_x  +  oId_patch_y  *  ctdb_in->patches_wide; 
int32  next,  fclass; 

feature  =  ctdb_intem_lookup_patch(ctdb_in,  patch_number); 
if  (feature) 

{ 

last  =  feature  +  feature->hdr.size; 
feature++; 

for(;  feature  <  last;  feature  +=  next) 

{ 

n_verts  =  feature->info.vertices; 
fclass  =  feature->info.feature_class; 
if  (fclass  ==  FC_MICRO_CMP) 

{next  =  NEXT_MICRO_CMP(n_verts);continue; } 
else  if  (fclass  =  FC_VOLUME_CMP) 

{next  =  NEXT_VOLUME_CMP(n_verts);continue; } 
else  if  (fclass  ==  FC_CANGPY_CMP) 

{next  =  NEXT_CANOPY_CMP(n_verts);continue; } 
else  if  (fclass  =  FC_AGGREGATE_CMP) 

{next  =  NEXT_AGGREGATE_CMP(n_verts);continue;  } 
else  if  (fclass  =  FC_LINEAR_CMP) 

{next  =  NEXT_LINEAR_CMP(n_verts);continue; } 
else  if  (fclass  ==  FC_LAID_LINEAR_CMP) 


This  section  skips  all  features 
except  the  roads,  rivers,  and 
other  linear  features. 
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next  =  NEXT_LAID_LINEAR_CMP(n_verts); 
if  (n_verts) 

{ 

if  (n_verts  >=  MAX_LAID_LINEAR_VERTS) 

{ 

printf("MAX_LAID_LINEAR_VERTS  (ct_cmplr.h)  exceeded.  Increase  it  to  %d\n", 
n_verts+l); 

n_verts  =  MAX_LAID_LINEAR_VERTS-1 ; 

} 

recompile_linears++; 
line.n_verts  =  n_verts; 
line. width  =  ctdb_in->incrD  * 

XYO_TO_POST(feature[2].laid_linear.  width, 0, 
ctdb_in->patch_incrD); 
line.soil  =  feature[l].ll_soiI.poly_char_index; 
if  (line.soil  ==  1 )  line. width  =  10.0; 
else  line. width  =  20.0; 
for  (i=0;i<n_verts;i++) 

{ 

xx  =  ctdb_in->incrD  * 

X  Y  0_T0_P0ST  (featu  re[3+i]  .xy.x,old_patch_x, 
ctdb_in->patch_incrD)  +  offset_x_meters; 
yy  =  ctdb_in->incrD  * 

XYO_TO_POST(feature[3+i].xy.y,old_patch _y, 
ctdb_in->patch_incrD)  +  offset_y_meters; 
line.verts[i][X]  =  xx; 
line.verts[i][Y]  =yy; 

} 

ctdb_buffer_laid_linear(&Iine); 

} 

} 

} 

} 

A.14  Conopy_triangles 

static  void  canopy_triangle  (micro, xO,yO,xl,yl,x2,y2) 
float64  xO,yO,xl,yl,x2,y2; 

CTDB_MC_POLY  GON  micro; 

{ 

micro.n_verts  =  3; 
micro.soil  =  720907; 
micro.reserved  =  0; 
micro.  verts[0][X]  =  x0; 
micro.verts[0][Y]  =  yO; 
micro.  verts[0][Z]  =  20.0; 
micro.verts[l][X]  =xl 
micro.verts[l][Y]  =yl; 
micro.verts[l][Z]  =  20.0; 
micro.verts[2][X]  =  x2 
micro.  verts[2][Y]  =  y2; 
micro.verts[2][Z]  =  20.0; 

} 

A.15  Collect  canopies 

static  void  collect_canopies  ( new_patch_x,new_patch_y,patch_x,patch_y, 


Roads  and  rivers  are  line  segments 
having  a  soil  type  and  a  width. 
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n_cpoIy,  canopy_polys,  n_cedge,  canopy_edges,cp_limit,  ce_limit) 
int32  patch_x,patch_y,new_patch_x,  new_patch_y; 
int32  *n_cpoly; 

CTDB_MC_POLY  GON  canopy_polysO; 
int32  *n_cedge; 

CTDB_CAN_EDGE  canopy_edges[|; 
int32  cpjimit,  ce_limit; 

{ 

int32  patch_number,  reference; 

CTDB_FEATURE_DATA_CMP  *feature,  *last,  *next; 

CTDB_MC_POLY GON  micro; 

CTDBCANJEDGE  edge; 
int32  verts,  fclass; 

int  ix,iy,ix_start,ix_stop,iy_start,iy_stop,my_flag; 
float64  x  1  ,y  1  ,x2,y2,d  1  ,d2,d3 ; 
float64  min_x,min_y,max_x,max_y; 
int  zero, one, two; 

patch_number  =  patch_x  +  patch_y  *  ctdb_in->patches_wide; 
feature  =  ctdb_intem_lookup_patch(ctdb_in,  patch_number); 
min_x  =  new_patch_x*ctdb_out.patch_incr*ctdb_out.incr; 
min_y  =  new__patch_y*ctdb_out.patch_incr*ctdb_out.incr; 
max_x  =  (new_patch_x  +  l)*ctdb_out.patch_incr*ctdb_out.incr; 
max_y  =  (new__patch_y  +  l)*ctdb_out.patch_incr*ctdb_out.incr; 
ix_start  =  min_x/mat_step;  iy_start  =  min_y/mat_step; 
ix_stop  =  max_x/mat_step;  iy_stop  =  max_y/mat_step; 
zero  =  0;  one  =  0;  two  =  0; 
verts  =  3; 

fclass  =  F  C_C  AN  OP  Y_CMP ; 
for  (iy  =  iy_start;  iy  <  iy_stop;  iy++) 

{ 

for  (ix  =  ix_start;  ix  <  ix_stop;  ix++) 

{ 

if  (canopy  [ix][iy]  ==  0 )  zero  =  zero  +  1; 
if  (canopy  [ix][iy]  ==  1  )  one  =  one  +  1; 
if  (canopy  [ix][iy]  —  2  )  two  =  two  +  1; 

} 

} 

if  ( (zero  —  0)  &&  (one  ==  0)) 

{ 

canopy J:riangIe(&micro,min_x,inin_y,max_x,inin_y,inax_x,inax_y); 
ctdbJbuffer_canopy(&micro); 

canopy_triangle(&micro,min_x,min_y,min_x,max_y,max_x,max_y); 

ctdb_buffer_canopy(&micro); 

} 

if  ( (two  >  0)  &&  (one  >  0)) 

{ 

reference  =  0; 

edge. fullness  =  DEFAULTFULLNESS; 
edge.start[X]  =  min_x;  edge.start[Y]  =  min_j; 
edge.end[X]  =  max_x;  edge.end[Y]  =  max_y; 
ctdb_buffer_cedge(&edge); 

if  ( canopy [ix_start]  [ iy_stop- 1  ]  ==  0  &&  canopy [ix_start]  [iy_start]  =  0 
&&  canopy[ix_stop-l][iy_start]  —  0  &&  canopy[ix_stop-l][iy_stop-l]  =  0  ) 

{ 

if  (  canopy [(ix_stop-l  +  ix_start)/2]  [iy_start  ]  —  2 
&&  canopy[(ix_stop-l+  ix_start)/2][iy_stop-l]  ==  2 
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&&  canopy[ix_start][(iy_stop-l  +  iy_start)  12]  =  0 
&&  canopy[ix_stop-l][  (iy_stop-l  +  iy_start)/2]  ==0) 

{ 

canopy_triangIe(&micro, 

(maxx  +  min_x)  /  4  +  (minjx  /  2),min_y, 

(maxx  +  min_x)  /  4  +  (minjx  /  2),min_y, 

(maxx  +  min_x)  14  +  (min_x  /  2),max_y); 
ctdb_buffer_canopy(&micro); 
canopy__triangle(&micro, 

(maxx  +  min  jx)  14  +  (minx  /  2),max_y, 

(maxx  +  min_x)  14  +  (min_x  /  2),min_y, 

(maxx  +  min_x)  14  +  (min_x  /  2),max_y); 
ctdb_buffer_canopy(&micro); 

} 

else 

{ 

canopy_triangle(&micro, 

(maxx  +  min_x)  /  2,  (max_y  +  min_y)  /  2, 
min_x,  (maxjy  +  min_y)  /  2, 

(maxx  +  minx)  /  2,max_y); 
ctdb_buffer_canopy(&micro); 
canopy_triangle(&micro, 

(maxx  +  minjt)  /  2,  (max_y  +  min_y)  /  2, 

(max  x  +  min_x)  /  2,  max_y 
max_x,(max_y  +  min_y)  /  2); 
ctdb_buffer_canopy(&micro); 
canopy_triangIe(&micro, 

(maxx  +  min_x)  /  2,  (max_y  +  min_y)  /  2, 

(max  x  +  minjx)  /  2,  min_y 
min_x,(max_y  +  min_y)  /  2); 
ctdbJbuffer_canopy(&micro); 
canopy_triangle(&micro, 

(max  x  +  min_x)  /  2,  (max_y  +  min_y)  /  2, 
max_x,(max_y  +  min_y)  /  2 
(maxx  +  min_x)  /  2,min_y); 
ctdb_buffer_canopy(&micro); 

> 

> 

else  if  ( canopy[ix_start][iy_stop-l]  —  0  &&  canopy[ix_stop-l][iy_start]  ==  2 ) 

{ 

reference  =  0; 

edge.fullness  =  DEFAULT_FULLNESS; 
edge.startfX]  =  max_x;  edge.start[Y]  =  min_y; 
edge.endfX]  =  max_x;  edge.end[Y]  =  max_y; 
ctdb_buffer_cedge(&edge); 

canopy_triangle(&micro,  min_x,  min_y,  max_x,min_y,  max_x,max_y); 

ctdb_buffer_canopy(&micro); 
if  ( canopy[ix_start][(iy_stop-l  +  iy_start)  /  2]  =  2  ) 

{ 

canopy_triangle(&micro,  min  x,  min_y, 
min_x,  (max_y  +  min_y)  /  2, 

(max_x  +  min_x)  /  2,  (max_y  +  min_y)  /  2); 
ctdb_buffer_canopy(&micro); 

} 

if  (  canopy[  (ix_stop-l  +  ix_start)/2  +  1][  (iy_stop-l  +  iy_start)/2  +2]  =  2  ) 

{ 
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canopy_triangle(&micro, 

(maxx  +  min  _x)  /  2,  max_y, 

(max_x  +  min_x)  /  2,  (max_y  +  min_y)  /  2, 

((max_x  +  min_x)  /  4)  +  (max_x  /  2,((max_y  +  min_y)  /  4)  +  (max_y  /  2)); 
ctdb_buffer_canopy(&micro); 

} 

if  ( canopy [  ((ix_stop-l  +  ix_start)  /  2)  -  1][  (iy__stop-l  +  iy_start)/2  +2]  ==  2  ) 

{ 

canopy_triang!e(&micro, 

(max_x  +  min  _x)  /  2,  (max_y  +  min_y)  /  2, 

((max_x  +  min  x)  / 4)  +  (min_x  /  2,  ((max_y  +  min_y)  / 4)  +  (max_y  / 2), 
(max_x  +  min_x)  / 2,  max_y); 
ctdbJbuffer_canopy(&micro); 

} 

} 

else  if ( canopy [ix_start][iy_start]  ==  0  &&canopy[ix_stop-l][iy_stop-l]  =  2) 

{ 

reference  =  0; 

edge.fullness  =  DEFAULT_FULLNESS; 
edge.start[X]  =  min_x;  edge.start[Y]  =  min_y; 
edge.end[X]  =  max_x;  edge.end[Y]  =  max_y; 
ctdb_buffer_cedge(&edge); 

canopy_triangle(&micro,  maxx,  min_y,  min  x,  max_y,  max_x,  max_y); 
ctdb_buffer_canopy(&micro); 

if  (  canopy[ix_start][  (iy_stop-l+  iy_start)/2  +  1]  ==  2 ) 

{ 

canopy_triangIe(&micro,  min_x,  (max_y  +  min_y)  /  2 , 

(maxx  +  minx)  /  2,  (max_y  +  min_y)  /  2, 
min_x,  max_y); 
ctdb_buffer_canopy(&micro); 

} 

if  ( canopy[ix_stop-l][iy_start+2]  ==  2  ) 

{ 

canopy_triangle(&micro,  (max_x  +  min_x)  /  2,min_y; 

(max_x  +  min_x)  /  2,(max_y  +  min_y)  /  2 
maxx,  min_y); 
ctdb_buffer_canopy(&  m  icro); 

} 

} 

else  if  (  canopy[ix_stop-l][iy_start]  ==  0  &&canopy[ix_start][iy__stop-l]  ==  2  ) 

{ 

reference  =  0; 

edge.fullness  =  DEFAULT_FULLNESS; 
edge.start[X]  =  min_x;  edge.start[Y]  =  min_y; 
edge.end[X]  =  min_x;  edge.end[Y]  =  max_y; 
ctdb_buffer_cedge(&edge); 
canopy_triangle(&micro,  min_x,  min_y, 
min  x,  max_y,  max_xmax_y); 
ctdbJbuffer_canopy(&micro); 

if  ( canopy[  (ix_stop-l  +  ix_start)/2  + 1][  (iy_stop-l  +  iy_start)  /  2]  ==  2  ) 

{ 

canopy_triangIe(&micro,  (max  x  +  min_x)  /  2,  (max_y  +  min_y)  /  2, 
max  x,  (max_y  +  min_y)  /  2,  maxjK,  max_y); 
ctdb_buffer_canopy(&micro); 

} 
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} 

else  if  (  canopy[ix_stop-l][iy_stop-l]  ==  0  &&  canopy [ix_start][iy_start]  —  2  ) 

{ 

reference  =  0; 

edge.fullness  =  DEFAULT_FULLNESS; 
edge.start[X]  =  min_x;  edge.start[Y]  =  min_y; 
edge.end[X]  =  min_x;  edge.end[Y]  =  max_y; 
ctdb_buffer_cedge(&edge); 

canopy_triangle(&micro,  min_x,  min_y,  min_x,  max_y,  max_x,  min_y); 

ctdb_buffer_canopy(&micro); 

if  (  canopy[ix_start+l][iy_stop-l]  =  2  ) 

{ 

canopy_triangle(&micro,  min_x,  max_y,  (max_x  +  min_x)  /  2,  max_y, 
(maxx  +  min_x)  /  2,  (max_y  +  min_y)  /  2); 
ctdb_buffer_canopy(&micro); 

} 

if  (  canopy[ix_stop-l][iy_start  +  2]  ==  2 ) 

{ 

canopy_triangIe(&micro,  (max_x  +  min_x)  /  2, 

(max_y  +  min_jy)  /  2,  maxx,  (max_y  +  min  _y)  /  2 
maxx,  min_y); 
ctdb_buffer_canopy(&micro); 

} 

} 

else  if  (  canopy [ix_start][iy_stop  - 1]  —  0 ) 

{ 

edge.start[X]  =  max_x;  edge.start[Y]  =  min_y; 
edge.end[X]  =  max_x;  edge.end[Y]  =  max_y; 
ctdb_buffer_cedge(&edge); 
canopy_triangle(&micro,  min_x,  min_y, 
max  x,  min_y,  max  x,  max_y); 
ctdb_buffer_canopy(&micro); 

} 

else 

{ 

edge.start[X]  =  max_x;  edge.startfY]  =  min__y; 
edge.end[X]  =  max_x;  edge.end[Y]  =  max_y; 
ctdb_buffer_cedge(&edge); 

canopy_triangle(&micro,  min_x,  min_y,  min_x,  max_y,  max_x,max_y); 
ctdb_buffer_canopy(&micro); 

} 


} 


A.16  Fit_vrt 

static  void  fit_vrt() 

{ 

int  khills,ix,iy,iix,iiy,icx,icy,number_grid_squares; 
float64  my_resi[500][500],x_center,y_center; 
float64  max_z,dx,dy,dxl,dyl,deg,my_z,loc_x,loc _y; 
float64  current_max,current_min; 
int  min_x,min_y,max_x,max_y; 
float64  len,len2,abs_z,rwx,rwy; 
int  isplit, iterations; 
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int  npx, npy, xstart, xstop, ystart, ystop; 

int  k,number_of_fits; 

max_z  =  -999999.0; 

npx  =  ctdb_in->max_x_post; 

npy  =  ctdb_in->max__y__post; 

for  (iy  =  0;  iy  <  npy;  iy++) 

{ 

for  (ix  =  0;  ix  <  npx;  ix++) 

{ 

my_resi[ix][iy]  =  original[ix][iy]  -  ctdb_in->min_z; 
temp[ix][iy]  =  0.0; 

abs_z  =  sqrt(my_resi[ix][iy]*my_resi[ix][iy]); 
if  ( abs_z  >  max_z) 

{ 

max_z  =  abs_z; 
icx  =  ix; 
icy  =  iy; 

} 

} 

} 

khills  =  0; 
isplit  =  2; 
while  (isplit  <  32) 

{ 

if  (isplit  <  8)  number_of_fits  =  4; 
if  (isplit  >=  8)  number_of_fits  =2; 
number_grid_squares  =  npx/isplit; 

x_center  =  (float)number_grid_squares*ctdb_in->incrD/2.0; 
for  (iiy  =  0;  iiy  <  isplit;  iiy++) 

{ 

ystart  =  iiy*number_grid_squares; 
ystop  =  (iiy  +  l)*number_grid_squares; 
if  (ystop  >  npy)  ystop  =  npy; 
for  (iix  =  0;  iix  <  isplit;  iix++) 

{ 

xstart  =  iix*number_grid_squares; 
xstop  =  (iix  +  l)*number_grid_squares; 
if  (xstop  >  npx)  xstop  =  npx; 
for  (k  =  0;  k  <  number_of_fits;k++) 

{ 

loc_y  =  (float)iiy*2.0*x_center  +  x_center; 
loc_x  =  (float)iix*2.0*x_center  +  x_center; 
hill[khills].hgt  =  0.0; 
hill[khills].wgty  =  0.4*drand48()  +  0.1; 
current_min  =  9999999.0; 
current_max  =  -9999999.0; 
for  (iy  =  ystart;  iy  <  ystop;iy++) 

{ 

for  (ix  =  xstart;  ix  <  xstop;ix++) 

{ 

if  (my_resi[ix][iy]  <  current_min) 

{ 

currentjnin  =  my_resi[ix][iy]; 
min_x  =  ix; 
min_y  =  iy; 

> 
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if  (my_resi[ix][iy]  >  current_max) 

{ 

current_max  =  my_resi[ix][iy]; 
max_x  =  ix; 
max_y  =  iy; 

} 

} 

} 

if  (  current_max*current_max  >  current_min*current_min) 

{ 

hill[khills].hgt  =  my_resi[max_x][max_y]; 
hill[khills].cx  =  (float)max_x*ctdb_in->incrD; 
hill[khills].cy  =  (float)max_y*ctdb_in->incrD; 
my_z  =  current_min; 

} 

else 

{ 

hill[khills].hgt  =  my_resi[min_x][min_y]; 
hill[khills].cx  =  (float)min_x*ctdb_in->incrD; 
hill[khills].cy  =  (float)min_y*ctdb_in->incrD; 
my_z  =  current_max; 

} 

hill[khiIls].power  =  (0.3  +  0.5*drand48()); 

Ien  =  sqrt((float)((max_x  -  min_x)*(max_x  -  min_x)  + 

(max_y  -  min_y)*(max__y  -  min_y))); 
len  =  len*ctdb_in->incrD; 
hill[khills].wgtx  = 

((xstop  -  xstart)*ctdb_in->incrD)*(0.7  +  0.3*drand48()); 
dx  =  (float)  (max_x  -  min_x); 
dy  =  (float)  (max_y  -  min_y); 
if  (dx*dx  >  0.001) 

{ 

hill[khills] .angle  =  90.0  -  atan(dy/dx)*45.0/atan(l  .0); 

} 

else 

{ 

hill[khills].angle  =  0.0; 

} 

hill[khills].hgt  =  hill[khills].hgt  +  ctdb_in->min_z; 
for  (iy  =  0;  iy  <  npy;  iy++) 

{ 

loc_y  =  iy*ctdb_in->incrD; 
for  (ix  =  0;  ix  <  npx;  ix++) 

{ 

loc_x  =  ix*ctdb_in->incrD; 

temp[ix][iy]  =  temp[ix][iy]  +  single_hill(loc_x,loc_y,khills); 
my_z  =  temp[ix][iy]  +  ctdb_in->min_z; 
my_resi[ix]fiy]  =  original  [ix]  [iy]  -my_z; 

} 

} 

khills  =  khills +1; 

} 

} 

} 

isplit  =  isplit*2; 

} 
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for  (nhills  =  khills;nhills  <  1000;nhills++) 

{ 

loc_x  =  icx*ctdb_in->incrD; 
loc_y  =  icy*ctdb_in->incrD; 
hill[nhills].cx  =  loc_x; 
hill[nhills].cy  =  loc_y; 

hill[nhills].hgt  =  my_resi[icx][icy]  +  ctdb_in->min_z; 

printfC'hill  #  %d  cx  %d  cy  %d  %.31f  \nM,nhiIIs,icx>icy,my_resi[icx][icy]); 

hill[nhills]. angle  =  180*drand48(); 

hillfnhills]  .power  =  0.7*drand48()  +  0.1; 

hill[nhills].wgtx  =  0. 1 0*drand48()*(ctdb _in->max_x  -  ctdb_in->min_x) 

+  0.005*(ctdb_in->max_x  -  ctdb_in->min_x); 

hill[nhills].wgty  =  0.8*drand48()  +  0.05; 

max_z  =  -999999.0; 

for  (iy  =  0;  iy  <  npy;  iy++) 

{ 

loc_y  =  iy*ctdbJn->incrD; 
for  (ix  =  0;  ix  <  npx;  ix++) 

{ 

loc_x  =  ix*ctdb_in->incrD; 

temp[ix][iy]  =  temp[ix][iy]  +  singIe_hiII(loc_x,loc_y,nhills); 
my_z  =  temp[ix][iy]  +  ctdb_in->min_z; 
my_resi[ix][iy]  =  original  [ix]  [iy]  -  my_z; 
abs_z  =  sqrt(my_resi[ix][iy]*my_resi[ix][iy]); 
if  (  abs_z  >  max_z) 

{ 

max_z  =  abs_z; 
icx  =  ix; 
icy  =  iy; 

} 

} 

} 

} 

} 

A.17  VRT 

static  double  vrt(px,py, khills) 
int  khills; 
double  px,py; 

{ 

double  my_z,current_term; 
int  ih; 

my_z  =  0.0; 

current_term  =  0.0; 

for  (ih  =  0;  ih  <=  khills;  ih++) 

{ 

currentterm  =  single_hill(px,py,ih); 
my_z  =  my_z  +  current_term; 

} 

my_z  =  my_z  +  ctdb__in->min_z; 
return  my_z; 

} 

A.18  Single_hiII 

static  double  single_hill(px,py,ih) 
int  ih; 
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double  px,py; 

{ 


/*  The  function  single_hill  evaluates  the  expression: 

h(px,py)=/zgre-w^Kx(Px’Py)-x(cx’cy))2+w^Y(px’Py)_Y^cx’cy^2^0Vi;er 

X(x,  y)  =  cos(co)x  +  sin(co)y  md  Y (x,  y)  =  -  sin(co)x  +  cos(co)y  v 

double  my_x,my_y,my_single_hill_z; 

double  xxx,yyy,xxjy}maxh,minh,scale_x,scale_y; 

double  cXjCyjCz^^s^O^OjdxXjdyyjexpon^z; 

double  aa,bb,my_w,my_sign?computed_wgt; 

double  current_term,offset_x,offset_y,hgt,hgtm; 

double  ndJength_of_axis; 

double  deg_toj*ad  =  atan(1.0)*4.0/180.0; 

double  my_max_z,my_m  in_z; 

int32  my_max_x?  myjmax_y,  my_min_x,  my_min_y; 

my_x  =  (px  -  ctdb_in->min_x)/(ctdb_in->max_x  -  ctdb_in->min_x); 

my_y  =  (py  -  ctdb_in->min_y)/(ctdb_in->max_y  -  ctdb_in->min_y); 

c  =  cos(deg_to_rad*hill[ih].angle); 

s  =  sin(deg_to_rad*hill[ih].angle); 

cx  =  (hill[ih].cx  -  ctdb_in->min__x)/(ctdb_in->max_x  -  ctdb_in->min_x); 
cy  =  (hill[ih].cy  -  ctdb_in->min_y)/(ctdb_in->max_y  -  ctdb_in->min__y); 
cz  =  (hill[ih].hgt  -  ctdb_in->min_z)/hgt_max; 
ndJength_of_axis  = 

(hill[ih].wgtx  -  ctdb_in->min_x)/(ctdb_in->max_x  -  ctdbjn->min_x); 
yyO  =  cx*c  +  cy*s; 
xxO  =  -cx*s  +  cy*c; 
my_sign  =  1.0; 

if  ( hill[ih].hgt  <  0.0)  my_sign  =  -1.0; 
aa  =  log(0.001/sqrt(cz*cz)); 

bb  =  log(ndjength_of_axis*ndjengthjrf_axis)*hill[ih].power; 

computed_wgt  =  -1.0*aa/exp(bb); 
yy  =  my__x*c  +  my_y*s; 
xx  -  -my_x*s  +  my_y*c; 
dxx  =  (xx  -  xx0)*(xx  -  xxO); 
dyy  =  (yy-yyO)*(yy-yyO); 

expon  =  exp(hilI[ih].power*log(dxx  +  hill[ih].wgty*dyy)); 
current_term  -  sqrt(cz*cz)*exp(-1.0*computed_wgt*expon); 
current_term  =  current_term  * hgt_max  ; 
my__single_hill_z  =  my_sign*current_term; 
return  my_single_hill_z; 

} 

A.  19  AddjSubcanopies 

static  void  add_subcanopies  (patch_water_state) 
int32  * patch_water_state ; 

{ 

CTDB_SEARCH_SPACE_PTR_CMP  search_space; 
int32  i,  j,ii,ij; 
int32  count; 
int32  my_type; 

CTDB_ABSTRACT_DATA_CMP  *data; 
int32  data_size; 


,  where 
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float32  *my_vert=  NULL; 
int32  vertexjiumber  =  0; 
float64  meters_perjpatch; 

meters_per_patch  =  ctdb_out.incrD  *  ctdb_out.patch_incrD; 
for  ( ii  =  0;ii  <  soil_region;ii-H-) 

{ 

vertex_number  =  my_reg[ii]  .count; 
f°r  ( jj  =  0;  jj  <  vertex_number;  jj++) 

{ 

my_vert[2*jj]  =  my_reg[ii].x[jj]; 
my_vert[2*jj  +  1]  =  my_reg[ii].y[jj]; 

} 

myjype  -  CTDB_ABSTRACT_SOIL_DEFRAG_CMP; 
data_size  =  1 ; 

data->soil_defrag.soil_index  =  9; 
data->soil_deffag.level  =  0; 
data_size  =  1; 

ctdb_store_abstract(&ctdb_out,  &quadjroot,  my_type, 
vertex_number,  my_vert,  data_size,  (int32  *)data, 
min_patch_x  *  meters_perj)atch, 
min_patch_y  *  meters_perjpatch, 
maxjatch_x  *  meters _per_patch, 
max_patch_y  *  meters_per_patch); 

} 

if  (numjverts  >  0)  free(lverts); 

} 
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A  SAMPLE  ENTITY  PARAMETER  FILE 
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A  SAMPLE  ENTITY  PARAMETER  FILE 


This  appendix  contains  the  entity  parameter  file  for  an  experimental  UGV.  The  file  resides 
in  the  directory  “/modsaD.O/common/src/ModSAF/entities,”  where  “/modsaf3.0/”  is  the  root 
directory  for  the  ModSAF  simulation  code.  In  the  following  text,  a  double  semi-colon  (;;) 
indicates  a  comment  in  the  file.  Some  formatting  and  bold  print  have  been  added  to  the  file  for 
clarity.  The  reader  is  referred  to  the  ModSAF  references  (Sagacitech  1997,  and  Smith  1995)  for 
more  information  about  entity  parameter  files. 


;  ;US_UGV_M_T_A _params.rdr 

US_UGV_M_T_A_MODEL_PARAMETERS  { 

;;  1.  ENTITY  PARAMENTERS 

(SMJEntity  DEFAULTJDEAD_RECKONING_PARAMETERS 

(vehicle_class  vehicleClassSimple) 

(guises  vehicle_US_UGV  vehicle_US_HMMWV) 
(send_dis_deactivate  true)  ) 

;;  2.  VULNERABILITY  MODELS  &  DAMAGE  ASSESMENT 
(filename  "dfdam_TRUCK.rdr") 

(damage_threshold  10.0)) 

(name  apcl)) 

(background  on) 

(sensors  commander-sight) 

(weapons ) 

VASSESS_ADA_GROUND 
VASSESS_ADA_THREATS 
VASSESS_IFV_OPTIONAL 
(no_target_load)  ) 

;;  3.  VEHICLE  COMPONENTS 

(SM_Components  (hull  SM_TrackedHull  SAFCapabilityMobility) 

(primary-turret  [SM_GenericTurret  |  0] ) 
(commander-sight  [SM_Visual  |  1] ) 

(gunner-sight  [SM_Visual  |  0] ) 

(radioA  [SM_GenRadio  |  0] ) 

(radioB  [SM_GenRadio  1 2] ) ) 

(SM_EnvAssess  (commander_sight  "commander-sight")) 

(SM_EnvReason ) 

;;4.  PATH  PLANNING  PARAMETERS 

(SMLocalMap  (skirt_deviation  0.3)) 

(SM_VMove  (background  on) 

(stopped_time  60.0) 

(default_speed  10.0) 

(defauIt_max_deviation  1000.0) 
(default_catchup_speed  0.0) 

(default_brake_strength  1.0) 

(max_backup_distance  1 .0) 

(planning_horizon  60.0) 

(execution_horizon  1 .0) 


(SM_DFDamage 

(SM_IFDamage 
(SM_V  Assess 
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(moving_obstacle_horizon  10.0) 
(env  sampling  period  6) ) 


;;  5.  TERRAIN  ANALYSIS  PARAMETERS 
(SM_Vterrain  (entity_period  100) 

(avoidance_mask  [VTERRAIN_BUILDING  | 
VTERRAIN_  WATER  | 
VTERRAIN_DITCH]  ) 


(avoid joils  SOIL_DEEP_WATER  123) 

(background  on) 

(movementjhreshold  1.0) 

(map_radius  500.0) 

(entity_radius  10.0) 

(history_list_spacing  20.0) 

(num_historyJist joints  50) 

(breach_obst_nominal_size  400.0)) 

;;  6.  MOBILITY  PARAMETERS 
(SM_TrackedHull  (mobility_model  1) 

(soils 

(SOIL_DEFAULT  (max_speeds  52.0  52.0)  SOIL_DEFAULT_TRACKED) 

(1 

(max_speeds  52.0  52.0) 
(max_decel  5.37) 
(max_climb  36.0) 

(max_accel  1.32) 

(max_tum  47.0) 

(dust_speeds  100.0  100.0  100.0) ) 

(2 

(max_speeds  52.0  52.0) 
(maxjiecel  5.37) 
(maxjlimb  36.0) 

(max_accel  1 .32) 

(max_turn  47.0) 

(dust_speeds  100.0  100.0  100.0) ) 

(3 

(max_speeds  52.0  52.0) 
(max_decel  5.37) 

(max  jlimb  36.0) 

(max_accel  1.32) 

(max_tum  47.0) 

(dust_speeds  100.0  100.0  100.0) ) 

(4 

(max_speeds  52.0  52.0) 
(max_decel  5.37) 
(max_climb  36.0) 

(max_accel  1.32) 

(max_turn  47.0) 

(dust_speeds  100.0  100.0  100.0) ) 

(5 

(max_speeds  52.0  52.0) 
(max_decel  5.37) 
(max_climb  36.0) 

(max_accel  1.32) 

(max_turn  47.0) 

(dust_speeds  100.0  100.0  100.0) ) 

(6 

(max_speeds  52.0  52.0) 
(max_decel  5.37) 
(max_climb  36.0) 

(max_accel  1.32) 

(max  turn  47.0) 

(dust__speeds  100.0  100.0  100.0) ) 

(7 

(max_speeds  52.0  52.0) 
(max_decel  5.37) 
(max_climb  36.0) 

(max_accel  1.32) 

(max_tum  47.0) 

(dust_speeds  100.0  100.0  100.0) ) 

(8 

(max_speeds  52.0  52.0) 
(max_decel  5.37) 
(max_climb  36.0) 

(max_accel  1.32) 

(max_tum  47.0) 

(dust_speeds  100.0  100.0  100.0) ) 

(9 

(max_speeds  8.05  8.05) 
(max_decel  5.37) 

(max_accel  1.32) 

(max_turn  47.0) 
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(max_climb  36.0) 

( 10  (max_speeds  52.0  52.0) 
(max_decel  5.37) 
(max_climb  36.0) 

( 1 1  (max_speeds  52.0  52.0) 
(max_decel  5.37) 
(max_climb  36.0)  . 

( 12  (max_speeds  52.0  52.0) 
(max_decel  5.37) 
(max_climb  36.0) 

(13  (max^speeds  52.0  52.0) 
(max_decel  5.37) 
(max__climb  36.0) 

( 14  (max_speeds  52.0  52.0) 
(max_decel  5.37) 
(max_ciimb  36.0) 

(15  (max__speeds  52.0  52.0) 
(max_decel  5.37) 
(max_climb  36.0) 

( 16  (max_speeds  52.0  52.0) 
(max_decel  5.37) 
(max_climb  36.0) 


(dust_speeds  100.0  100.0  100.0) ) 

(max_accel  1 .32) 

(max_tum  47.0) 

(dust_speeds  100.0  100.0  100.0) ) 

(max_accel  1.32) 

(max_turn  47.0) 

(dust_speeds  100.0  100.0  100.0) ) 

(max_accel  1.32) 

(maxjurn  47.0) 

(dust_speeds  100.0  100.0  100.0)) 

(max_accel  1.32) 

(max_tum  47.0) 

(dust_speeds  100.0  100.0  100.0)) 

(max_accel  1.32) 

(maxjurn  47.0) 

(dust_speeds  100.0  100.0  100.0) ) 

(max_accel  1.32) 

(maxjurn  47.0) 

(dust_speeds  100.0  100.0  100.0)) 

(max_accel  1.32) 

(max_tum  47.0) 

(dust_speeds  100.0  100.0  100.0) ) ) 


(fueljisage  (0.0  100.0)  (0.125  12.0))) 

;;  7.  TURRET  PARAMETERS 

( [SM_GenericT urret  |  0]  (physdb_name  "primary-turret") 

(rates  continuous  0.0  40.0)) 

;;  8.  VISUAL  SYSTEMS 

([SM_Visual  |  0]  VISUAL>PC_DRIVER_DVO^NVO ) 

([SM_Visual  |  1]  VISUAL_LOSAT_HIRESJR ) 

(SM_SubComp) 

(SM_VSpotter  (background  on) 

(sensors  commander-sight) 

VSPOTTER_SPECS  ) 

;;  9.  RADIOS 

([SM_GenRadio  |  0]  (net_name  "platoon jiet") 

(aspid  ASPID_MODSAF_TEXT) 
GENRADIO_BLUE_P ARAMS  ) 

([SM_GenRadio  |  2]  (net_name  "company_net") 

(aspid  ASPID_MODSAF_TEXT) 

genradio_blue_params  ) 

;;  10.  VISUAL  SEARCH  ALGORITHMS 
(SMJVsearch  (search_type  ground) 

(scan_mode  static) 

(background  on) 
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(turret_scanner  "primary-turret"  9.0  5.0) 
(gun_scanner  "") 

(visual_scanners  "commander-sight") 
(stopped_duty_cycle  1.0) 
(moving_duty_cycle  0.0) 

(restrict2for  "none")) 


} 
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1992).  The  new,  higher  resolution  terrain  directly  affects  the  vehicle  dynamics  and  the  line-of-sight  algorithms.  The  new  terrain 
does  not  directly  affect  the  ModSAF  route-planning  algorithms.  In  the  second  phase  of  our  terrain  database  modifications,  the 
slopes  on  the  new  terrain  are  examined.  Regions  that  are  steep  or  inaccessible  to  the  XUV  are  marked  as  obstacles  in  the  database. 
The  route-planning  algorithms  use  these  “obstacles”  to  avoid  planning  routes  through  regions  that  are  too  steep  for  the  XUV. 
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